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ABSTRACT 


A  system  has  been  developed  using  the  Navy  Naviga¬ 
tion  Satellites  to  provide  azimuth  information  on  a  worldwide 
basis.  This  capability  is  particularly  useful  in  the  polar  re¬ 
gions  where  conventional  navigation  systems  based  on  gyro¬ 
compasses  are  handicapped  by  gyro  alignment  difficulties. 

In  this  newly  developed  system,  two  antennas  separated  by 
50  to  100  meters  a^e  used  as  an  interferometer  array  to 
measure  the  phase  difference  of  a  satellite  signal.  Based 
on  this  information,  it  is  possible  to  determine  the  azimuth 
of  the  baseline  connecting  the  two  receiving  antennas. 
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FOREWORD 


Precision  navigation  has  been  a  problem  ever  since 
man  began  to  travel  into  unfamiliar  areas.  Today  the  prob¬ 
lem  is  crucial  whether  it  is  related  to  driving  a  car  on  an 
unfamiliar  road  or  piloting  an  airplane  across  an  ocean, 
where  no  markers  exist. 

The  goal  of  navigation  can  be  simply  stated  as  suc¬ 
cessful  travel  from  an  arbitrary  location  to  a  desired  loca¬ 
tion  in  an  optimal  manner.  To  solve  this  problem  the  navi¬ 
gator  needs  two  pieces  of  information:  his  current  location 
and  a  reference  direction.  (It  is  assumed  that  the  navigator 
knows  the  location  of  his  destination  and  all  intervening  dis¬ 
turbances.  )  To  be  useful,  this  information  should  be  readily 
available  over  the  total  area  of  possible  travel. 

The  reference  direction  can  be  defined  in  terms  of 
azimuth.  Azimuth  is  an  angle  referenced  to  the  earth's 
spin  axis  such  that  an  angle  of  0°  is  called  north,  east  is 
90°,  south  is  180°,  and  west  is  2  70°. 

A  solution  to  the  problem  of  determining  a  reference 
direction  will  be  treated  here.  The  proposed  technique  uses 
the  radio  signal  received  from  a  satellite  whose  location  is 
known  as  a  function  of  time  a. id  tw  o  nondirectional  antennas 
used  as  an  interferometer.  This  system  has  been  demon¬ 
strated  using  t lie  Navy  navigation  satellites. 

The  interferometer  method  oilers  many  advantages 
over  other  "north  seeking"  techniques.  The  most  common 
of  these  methods  are  starsighting  and  north- seeking  gyro 
compasses.  Starsighting  systems  are  weather  limited  and 
primarily  useful  at  night.  N  orth- seeking  gyro  compasses, 
on  the  other  hand,  overcome  the  weather  and  night  limitations 
but  w  ill  not  self- align  in  the  high  latitude  polar  regions  and 
are  subject  to  drift  errors.  The  interferometer  azimuth  de¬ 
termination  system  overcomes  the  disadvantages  of  both  of 
these  methods  by  allowing  measurements  on  a  worldwide 
basis,  independent  of  weather  or  time  of  day. 
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1.  INTRODUCTION 


THE  NAVY  NAVIGATION  SATELLITE  SYSTEM 

The  Navy  Navigation  Satellite  System  (NAVSAT)  de¬ 
veloped  by  The  Applied  Physics  Laboratory  for  the  Depart¬ 
ment  of  the  Navy  consists  of  four  or  more  near-earth  satel¬ 
lites  in  polar  orbits,  tracking  stations,  injection  stations,  a 
computing  center,  and  receiving  navigation  equipment.  This 
system  employs  doppler  tracking  for  both  satellite  position 
determination  and  navigation.  To  determine  satellite  posi¬ 
tion,  four  tracking  stations  in  precisely  known  locations  ob¬ 
serve  t lie  doppler  shift  of  the  highly  stable  radio  signals  gen¬ 
erated  by  the  satellite  transmitter  as  the  satellite  passes  over 
the  station.  The  doppler  information  allows  the  computing 
center  to  determine  the  satellite's  position  as  a  function  of 
time.  From  this  measured  positional  information,  the  loca¬ 
tion  of  the  satellite  at  a  future  time  can  be  predicted.  These 
predictions  are  stored  in  the  memory  of  the  satellite  by  the 
injection  station.  As  the  satellite  orbits  t lie  earth,  it  contin¬ 
ually  transmits  its  position  and  precise  time.  By  receiving 
time  and  position  data  from  the  satellite  and  observing  the 
doppler  shift  in  the  satellite  signals  a  navigator  can  determine 
his  position. 

The  Navy  navigation  satellites  are  in  nominal  600  nmi 
circular  polar  orbits  with  periods  of  approximately  105  min¬ 
utes.  There  are  currently  five  such  satellites  with  approxi¬ 
mately  equal  angular  separations  between  the  orbital  planes 
(Fig.  1)  so  that  the  coverage  is  fairly  evenly  spread,  pro¬ 
viding  reception  possibilities  approximately  once  every  2 
hours  at  the  equator  and  considerably  more  frequently  as  lat¬ 
itude  incre'etses.  Typically  the  satellite  is  in  view  for  approx¬ 
imately  15  minutes  during  each  pass. 

The  ground  support  system  consists  of  tracking  sta¬ 
tions  that  receive,  record,  and  digitize  doppler  signals  from 
the  satellites,  a  computing  center  that  computes  future  orbits, 
orbital  parameters,  and  time  corrections;  and  an  injection 
station  that  transmits  these  new  orbital  parameters  and  time 
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NAVY  NAVIGATION  SATELLITE  SYSTEM 


Fig  1  THE  NAVY  NAVIGATION  SATELLITES  IN  POLAR  ORBIT 
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corrections  to  the  satellite.  In  addition,  the  U.S.  Naval  Ob¬ 
servatory  recovers  the  satellite  time  signals  and  compares 
them  with  universal  time.  This  information  is  forwarded  to 
the  computing  center  for  the  time  correction  computations. 

Navigation  satellite  data  can  be  used  to  aid  navigation 
systems  in  a  number  of  ways,  such  as  providing  survey  bench 
marks,  or  to  update  the  inertial  platforms  of  ships  or  air¬ 
craft.  These  navigation  tasks  can  be  performed  on  a  world¬ 
wide  basis  regardless  of  weather. 


THE  AZTRAN  SYSTEM 

Aztran,  an  abbreviation  for  "azimuth  from  Transit,  " 
is  a  system  that  determines  north  reference  from  a  single 
pass  of  a  navigation  satellite.  The  system  uses  two  antennas 
separated  by  a  baseline  distance  of  30  to  200  meters.  Be¬ 
cause  of  the  separation  of  the  antennas,  there  is  a  transmis¬ 
sion  path  length  difference  from  the  antennas  to  the  satellite. 
From  Fig.  2  it  can  be  seen  that  the  path  difference  m ay  be 
computed  from  the  equation: 

Path  difference  =  D  •  cos  (ET  )  •  cos  (Az  -Az  ). 

L  s  at  ant 

When  the  signals  from  the  two  antennas  are  compared,  there 
is  a  phase  difference  caused  by  this  path  length  difference: 

Phase  difference  =  (2TTD/X)  -  cos  (ET  )  cos  (Az  -Az  ), 

L  sat  ant 

where  X  is  the  wavelength  of  the  received  signal. 

As  the  satellite  moves  across  the  sky,  there  are 
changes  in  the  path  length  difference  and  the  corresponding 
phase  angle.  This  phase  change  is  measured  between  time 
A  and  time  B  and  is  equal  to: 


2  ff  d  .  r  , 

°  =  \ —  toos  (E  )  cos(Az  -Az  )j  -Icos(E  )cos(Az  -\z  )],} 

A  L  Sat  ant  B  L  sat  ant  ^ 
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In  this  equation  (since  the  observer's  and  satellite's  positions 
are  known)  the  only  unknowns  are  D,  the  antenna  separation, 
and  the  antenna  baseline  azimuth  (Azan^).  By  making  repeated 
measurements  over  fixed  time  intervals  during  a  satellite 
pass,  the  baseline  azimuth  and  antenna  separation  are  deter¬ 
mined. 


APPLICATIONS 

Possible  applications  of  such  a  system  are  numerous. 
Accurate  azimuth  information,  easily  obtained  in  the  field, 
could  be  used  for  gun  fire  control  and  for  mobile  missile 
field  sites.  This  capability  would  free  the  missile  crews 
from  their  dependence  on  presurveyed  sites  for  this  informa¬ 
tion  and  guarantee  its  availability  day  or  night,  regardless 
of  weather.  Another  possible  application  is  in  the  alignment 
of  inertial  platforms.  Aztran  would  be  particularly  useful 
in  the  polar  regions  where  an  inertial  platform  will  not  self- 
align  and  must  be  referenced  to  an  external  azimuth  source. 
Ships  at  sea  and  aircraft  are  other  potential  users. 
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2.  THEORY  OF  OPERATION 


ORBIT  CALCULATIONS 

The  Navy  navigation  satellites  transmit  a  message 
containing  the  predicted  Keplerian  elements  of  the  satellite's 
orbit  in  inertial  space,  and  a  precise  time  mark  every  even 
2  minutes.  These  Keplerian  elements  are  then  converted 
into  Cartesian  earth-centered  coordinates,  which  rotate  with 
the  earth  such  that  a  fixed  point  on  the  earth  has  fixed  x-y-z 
coordinates. 

Consider  the  ellipse  (PSA)  shown  in  Fig.  3  with  a  cir¬ 
cle  (PCA)  circumscribing  the  ellipse.  Point  O'  is  the  center 
of  the  circle  and  the  ellipse.  Point  O  is  at  one  focus  of  the 
ellipse  which  is  the  origin  of  the  u-v  coordinate  system  and 
the  geocenter.  The  distance  AQ  is  called  the  semi- major 
and  BQ  is  called  the  semi- minor  axis.  Also: 


where  €  is  called  the  eccentricity  of  the  ellipse  and  can  vary 
from  0  (for  a  circle)  to  1  (for  a  straight  line). 

If  this  ellipse  represents  the  orbit  in  the  u-v  plane 
then  O  is  the  center  of  the  earth  and  S  is  the  position  of  the 
satellite  at  some  arbitrary  time  t.  The  t'me  at  which  the  sat¬ 
ellite  is  at  its  lowest  point  or  perigee  (point  P)  is  called  tp. 

The  angle  E  as  shown  in  Fig.  3  is  called  the  accentric 
anomaly.  The  mean  (average)  motion  of  the  satellite  is  de¬ 
fined  as: 


2  IT 

n  =  -j-  . 

where  T  is  the  orbital  period  of  the  satellite.  The  mean  angu¬ 
lar  motion  (M)  is  the  average  angular  motion  in  radians  since 
the  last  time  of  perigee: 

M  =  —  •  (t  -  t  )  -  n  •  (t  -  t  ) 

T  p  p 
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Therefore,  given  AQ,  €,  n,  and  tp,  the  u- v  coordinates 
of  the  satellite  at  S  can  be  found  as  a  function  of  time  as  fol¬ 
lows- 


M  (t)  =  n  •  (t  -  t  ) 

P 

E  (t)  =  M  (t)  +  C  sin  M  (t)  . 

Then: 


u(t ) 

Aq  *  [cos  E(t)  -  f] 

v(t) 

A  •  V 1  -  «  2  sin  E  (t) 
o 

Equation  (1)  would  hold  true  for  an  unperturbed  orbit 
about  a  spherical  mass.  The  orbit  of  an  earth  satellite  does 
not  remain  in  an  unperturbed  orbit  due  to  variations  in  the 
gravitational  field,  air  drag  on  the  satellite,  and  the  effects 
of  other  solar  bodies  (primarily  the  sun  and  moon).  There¬ 
fore,  some  of  the  orbital  parameters  which  were  treated  as 
constants  must  be  made  functions  of  time,  and  cross-plane 
terms  must  be  added. 

Consider  making  the  u-v  plane  into  a  right-hand  Car¬ 
tesian  coordinate  system  by  the  addition  of  a  w  axis.  Out- 
of-plane  variations  are  then  solely  in  this  axis. 

In  the  Xavy  navigation  satellite  system  the  semi-major 

axis  and  the  eccentric  anomaly  are  treated  as  functions  of 

time.  A  cross-plane  term,  T?(t ),  is  also  introduced.  The 

coordinates  of  the  satellite  in  the  u-v-w  frame  can  now  be 

found  as  follows,  given  A  ,  LA( t),  €,  n,  t  ,  AE(t),  and  T](t): 

o  p 

M(t)  =  n(t  -  t  ) 

P 

E(t )  =  M(t)  +  €  sin  M(t )  +  iE(t) 

A(t )  =  A  +  ^A(t) 
o 
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"u(t)“ 

A ( t )  •  [cos  E(t)  -  cj 

v(t ) 

= 

A(t )  •  x/ 1  €  “  sin  E(t) 

_\V  (t)_ 

_Tj(t  ) 

The  u-v-w  coordinate  system  can  now  be  transformed 
into  an  earth- centered  non  rotating  coordinate  system.  De¬ 
fine  the  X  -V  -Z '  coordinate  system  as  having  the  origin  at 
the  center  of  t lie  eartli  and  the  positive  Z  axis  pointing  north 
along  the  earth's  spin  axis.  The  X/-Y/  plane  lies  in  the  equa¬ 
torial  plane  and  the  positive  X'  axis  points  toward  Aires  (y). 
This  coordinate  frame  is  fixed  in  inertial  space. 

Figure  4  shows  the  u-v-w  and  the  X'-Y'-Z'  coordi¬ 
nates.  The  origins  of  the  two  systems  coincide  and  the  v 
axis  lies  in  the  X/-Y  plane.  The  longitude  of  the  ascending 
node  (0 0)  is  the  angle  measured  in  the  equatorial  plane  from 
the  positive  X  axis  to  the  negative  v  a.X’3.  The  orbital  incli¬ 
nation  (i)  is  the  angle  measured  from  the  positive  Z'  axis  to 
the  positive  w  axis.  The  argument  of  perigee  (to)  is  the  angle 
from  the  ascending  node  to  t lie  point  of  orbital  perigee. 

Given  the  preceding  parameters  defining  the  satellite's 
orbit,  Eq.  (3)  is  the  transformation  from  the  u-v-w  coordi¬ 
nate  system  to  the  X'-Y'-Z7  system: 


\  ‘ 

rnsif!)  si'  (SI )  0 

10  0 

(  os(*' 

-  sin(u.‘)  0 

u 

'l 

strilfi  )  i  ns(D  0 

0  cost  i)  -  sin(i ) 

sin(  jj  ) 

<  os(u.’)  0 

V 

y 

0  0  1 

0  si’  (i )  r os( i ) 

0 

0  1. 

w_ 

Figure  5  shows  the  re lationship  between  the  X'-Y'-Z7 
coordinate  system  and  the  desired  X-Y-Z  frame  which  is  fixed 
with  respect  to  the  earth's  surface.  The  origins  and  Z  axis 
of  the  two  systems  are  coincident.  cc  is  the  earth's  angular 
rotation  rat”.  The  X-Y  plane  and  the  X/-Y/  plane  are  both 
equatorial  and  rotated  by  \g,  the  Greenwich  hour  angle.  If 
the  angle  /3  is  defined  as  Ag  -  cce  (t  -  tp),  then: 
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V 

cos(^)  -sin(0)  0 

x'" 

y' 

= 

-sin(/3)  cos(£)  0 

Y ' 

z 

0  0  1 

_Z'_ 

(4) 


SATELLITE  MESSAGE 


During  each  2  minute  interval  a  Navy  navigation  satel¬ 
lite  transmits  19  message  words  used  for  navigation.  The 
first  8  are  called  ephemeral  words  and  contain  variables  for 
the  3  past  2  minute  intervals,  the  current  interval,  and  the 
next  4  intervals.  The  variables  contained  in  each  ephemeral 
are: 

t.  -  the  time  of  the  interval, 
k 

jE^  -  the  variation  in  the  eccentric  anomaly. 


-A^  '  the  variation  in  the  semi-major  axis,  and 
-  the  out  of  plane  distance. 

The  next  11  words  received  do  not  change  with  each  2  minute 
interval.  They  are: 

t  -  time  of  satellite  perigee, 

n  -  mean  motion  of  the  satellite, 

u;  -  the  argument  of  perigee  at  t  , 

P 

x  -  precession  rate  of  perigee, 

€  -  orbital  eccentricity, 

Aq  -  semi-major  axis, 

H  -  right  ascension  of  ascending  node  at  t  , 

o  p 
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-  precession  rate  of  node, 

Ci  -  cosine  of  the  inclination, 

A  -  longitude  of  Greenwich  at  t  ,  and 

g  P 

Si  -  sine  of  the  inclination. 


Figure  6  is  an  example  of  a  sequence  of  2  minute  mes¬ 
sages  received  trom  the  satellite.  A  minimum  of  three  in¬ 
tervals  is  required  to  solve  the  navigation  equations  for  the 
latitude  and  longitude  of  the  observer  and  the  difference  in 
the  frequencies  of  the  satellite  and  receiver  oscillators. 

This  last  quantity  is  not  predictable  and  must  be  a  variable 
ot  the  solution.  A  satellite  pass  of  maximum  duration  will 
be  observable  for  seven  or  eight  2- minute  intervals. 


The  set  of  orbital  elements  is  then  arranged  as  fol¬ 
lows,  where  k  represents  a  2  minute  interval  during  the  pass, 
and  Tk  is  the  correct  time  at  a  2  minute  time  mark: 


t 


k 


M 


E 


k 


A 


k 


=  M,  +  *  sin  M,  + 
k  k  k 

=  A  +  hA 
o  k 


Uk 

\  ■ 

(cos  E.  - 
k 

«) 

Vk 

V 

Vi  -  *2 

sin  E. 
k 

_"k_ 

> 

Equation  (3)  is  then  transformed  into  the  X-Y-Z  frame  using 
Eq.  (3)  and  (4)  for  each  interval  where: 


a:.  -  a:  -  a'  •  t, 

k  o  k 

Ft  =  n  +  r?  •  t 

k  o  k 
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420050084 
23  03  I  I  €2 0 
800015910 
<0999991  0 
43  00  I02"9 
2492  9?  I  I  0 
«.33l  5440 
90999991  0 
640109470 
2532  52  40  4 
99  3  0  1  5  4  4  3 
90  9  9  9  9  9  1  0 
2301 93757 
2631  «2684 
<0901  5640 
"09999910 
210261073 
270052  995 
90301  5440 
•09999910 
220501413 
090323926 
930015640 
9099999 1  0 
23031  I  775 
090133075 
930015640 
809999910 


6300  70213 
240261  955 
907463490 

643100470 
253252  425 
907463490 

200190757 
2601  924.84 
907463480 

210261070 
2  70092495 
607463490 

220301 413 
0  9  002  3  026 
807465480 

23031 1775 
090133075 
807463460 

2402  92  110 
400243049 
807463480 


6401 70404 
029971340 
8343.1  8690 

200190757 
094036060 
8|  43  I  6990 

210261070 

094036060 

834316990 

220301 413 
094036060 
934316990 

23031  I  775 
094036060 
834316980 

2402  921  10 
094036060 
834316980 

250252425 
094036060 
83  4316  980 


200250662 
8366  13670 
900002660 

2  I  02  6 1 0  70 
836612930 
9000023  80 

22  0301  415 
83 6C  I  293  0 
90  0  0  02  3  8  0 

230311 775 
836612930 
900002380 

2402921  10 
836612930 
900002380 

250252425 

836612930 

900002380 

260182684 

836612930 

900002380 


210300960 

934694710 

800062340 

220301413 

835553830 

800062390 

23031  I  775 
835553930 
800062390 

2  402921  10 
835553930 
800062390 

250252425 

834553930 

800062390 

260182684 

834553830 

800062390 

2  70092895 
835553830 
800062390 


220321289 
8001 99950 
809826660 

23031 1775 
8001 99470 
82591  I  760 

240292  I  10 
8001 95470 
92  591  1760 

250252  425 
8001 99470 
82591 1760 

260182684 
800199470 
82591 1760 

2  70092895 
800199470 
82591 1760 

080023026 

800199470 

825911760 


VCTEO  MESSAGE 


2001  90757 
260182684 
800199470 
82591  I  760 


210261070 
2  70092895 
800015640 
809999810 


220301413 

080023026 

907463480 


23031  I  775 
094036060 
834316980 


2  402  921  10 
836612930 
900002380 


250252425 

835553830 

800062390 


000243  777433  307700 
000253  630523  654600 
000263  I  12623  250400 
000273  777777  777740 


000000  055300  541002 
751474  740100  051350 
7T7777  777740  777777 
777777  777740  777777 


213700  557250  611440 
117600  102243  105100 
777740  777777  "’77740 
777740 


DAY .  288 

TIME....  1430 


Fig.  6  EXAMPLE  OF  SEQUENCE  OF  2  MINUTE  SATELLITE  MESSAGES 
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More  information  concerning  the  details  of  the  solu¬ 
tion  of  the  navigation  equations  can  be  found  in  Technical 
Memorandum  TO  819-1,  Program  Requirements  for  Two- 
Minute  Integrated  Doppler  Satellite  \avigation  Solution  edited 
by  J.  B.  Moffett. 

AZIMUTH  EQUATION’S 

In  this  section  the  equations  are  developed  which  allow 
the  determination  of  the  azimuth  of  a  line  connecting  two  an¬ 
tennas.  The  relative  positions  of  the  two  antennas  can  be  in¬ 
terpreted  in  terms  of  azimuth  angle  or  in  earth-fixed  coordi¬ 
nates. 


Recall  the  earth  fixed  X-Y-Z  coordinate  system  from 
the  previous  section  where  the  center  of  the  earth  was  the  or¬ 
igin,  the  equator  lay  in  the  X-Y  plane  and  the  Greenwich  me¬ 
ridian  lay  in  the  X-Z  plane  (Fig.  7).  If  O  is  the  geodetic  lati¬ 
tude  and  X  is  the  geodetic  longitude,  then: 


X 


Y 


Z 


_  I  _ Re _ 

n  j  -  2  .,2  ,2  -  f 

I  ,  cos  o  +  ( 1-  f)  sin  n  j 

< _  Re 

n  '  -  2  ..2  .  2  -Y 

I  .  COS  o  +  (  1-  f)  Sill  Oj 

_  ( _ Re(l-f)2 _ 

n  9  9  2  i 

|  _cos“  O  +  ( 1  -  f)“  sin  o  j 


+  Ut 


+  Ht 


+  lit 


cos  o  cos  X 


cos  o  sin  X 


sin  O  , 


(6) 

(7) 

(8) 


where  (X11,  Yn,  Zn )  are  the  coordinates  of  the  navigator  on 
the  earth's  surface  whose  latitude  and  longitude  are  respec¬ 
tively  o  and  X,  Re  is  the  equatorial  radius  of  the  earth,  lit 
is  the  observer's  height  above  the  reference  ellipsoid,  and 
f  is  t he  coefficient  of  flattening  of  the  reference  ellipsoid. 
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Assume  that  the  position  of  a  satellite  is  known  in  the 
earth  fixed  coordinate  frame  and  is  (Xs,  Ys,  Zg).  Then  the 
slant  range  from  satellite  to  navigator  can  be  expressed  as: 


SR  =  [  (x  -X  )2  +  (Y  -Y  )2  +  (Z  -Z  )2 
s  n  s  n  s  n 


0  O  9  *. 

=  (x  “  +  y  “  +  z  “r  , 

net  net  net 


where 


(9) 


X 

net 


=  (X  -Y  ),  etc. 
s  n 


Let  the  second  antenna  be  located  at  (Xn  +  AX,  Y  +  AY,  Zn 
+  *Z)  (Fig.  8),  at  an  azimuth  angle  0  from  North,  and  with 
an  antenna  separation  D.  Then: 


2  o  2  * 

D  =  (  _X  +  AY"  +  AZ  )c  .  (10) 

The  difference  in  slant  range  from  the  satellite  to 
the'  two  antennas  (-SR)  can  then  be  calculated: 


SK  |  [Xs  -  (\(+  X)j2  +  [VH  -  (Y  +’.Y)]“  +  [Z_  -  (Z  +  Z)]2!  *  -  S 


s  n 


s  n 


SR 


not 


4-  [V 


not 


'A  ■ 


(11) 


Expanding  this  results  in: 


-SR  =  [SR"  -  2(X  AX  +  Y  AY  +  Z  AZ)  +  D2]^  -  SR  . 
net  net  net 


Neglecting  (D“)  gives: 


ASR  =  SR  [  1  -  (X  ^  AX  +  Y  AY  +  Z  aZ )f  -  SR  . 
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I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


Taking  the  first  term  of  the  binominal  expansion  gives  the 
following: 


ASR  =  SR  -  —  (X  AX  +  Y  AY  +  Z  AZ)  -  SR 
SR  net  net  net 


ASR  =  - 


1 

SR 


AX  + 


Y  AY  +  Z  AZ)  . 
net  net 


(12) 


The  only  remaining  problem  is  to  determine  AX,  AY, 
and  AZ  in  terms  of  the  antenna  separation  and  azimuth.  Let 
the  antenna  separation  be  D,  the  azimuth  angle  be  0,  and  the 
difference  in  antenna  heights  above  the  reference  ellipse  be 
AH 


The  in-plane  distance  between  the  antenna  locations 
^plane^  *s: 


D„la„e  *  <°2  - 


From  Figs.  9  and  10  and  the  preceding  definitions,  it 
can  be  shown  that: 


~x 

=  -I)  , 

(sin  o 

cos  X  c°s 

;  0  +  sin  X 

sin  9  )  + 

M  cos  o  cos  X 

(13) 

plane 

n 

n 

n  n 

n 

n  n 

-Y 

=  -D  . 

(sin  o 

sin  X  cos 

9  -  cos  X 

sin  9  )  + 

H  cos  sin  X 

(14) 

plane 

n 

n 

n  n 

n 

n  n 

-Z 

--  I)  .  c 

os  o  cos  6  +  fill 

sin  O 

(15) 

plane 

11 

n 

n 

These  expressions  for  AX,  AY,  and  AZ  can  be  substi¬ 
tuted  back  into  Eq.  (12),  giving  : 

i>  . 

SR  P, 7  r  u  (\  cos  X  +  V  sin  A  )  sin  O  -  X  cos  (P  j  cos  6 

SR  net  net  net 

+  (X  sin  A  -  V  cos  A  )  sin  6  ] 
net  net 

-  4r-  (X  cos  A  +  Y  sin  A )  cos 0  +  7.  sin  <P  ]  .  (16) 

SR  net  net  net 
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The  path  length  difference  between  the  two  antennas 
(ASR)  can  be  expressed  in  wavelengths  as  NX,  where  X  is  the 
wavelength  of  the  frequency  being  received  and  N  is  the  num¬ 
ber  of  cycles.  If  another  measurement  were  taken  at  a  later 
time,  ASR2  =  N2X. 

It  is  not  convenient  to  determine  the  actual  value  of 
N.  Therefore,  the  change  in  phase  with  time  is  determined 
by  measuring  the  phase  change  between  times  tj  and  t2- 
This  gives: 

uSR2  -  -SR  =  (N2  -  N  >X  . 

The  process  is  repeated  far  as  many  time  intervals  as  possi¬ 
ble  during  a  satellite  pass  in  a  manner  similar  to  the  satel¬ 
lite  location  technique.  Using  this  method,  it  never  becomes 
necessary  to  know  the  original  N^;  knowledge  is  required  of 
only  the  phase  differences. 


APPLICATIONS 


There  are  a  variety  of  solution  forms  to  the  problem, 
depending  on  what  information  is  known  and  what  is  needed, 
as  listed  below: 


Possible 

Application 


Known 

Variables 


Solution 

Variables 


Minimum 
Required 
Number  of 
Intervals 


Directions  in  earth- 
fixed  coordinates 

None 

AX,  AY,  AZ 

for  aligning 

D 

AX,  AZ  or 

navigation  platforms 

AX,  AY  or 

AY,  AZ 

Azimuth  determination 

'-‘H 

A  '.imuth  (6) 

ah,  d 

Azimuth  (9) 

3 

2 

2 

1 
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Because  the  quantities  being  measured  in  each  inter 
val  are  not  exact,  but  are  corrupted  by  noise,  it  is  desirable 
to  obtain  more  data  intervals  than  required.  It  will  be  shown 
in  a  later  section  how  dramatically  the  results  improve  w'ith 
increasing  number  of  data  points. 

Assume  that  one  wishes  to  solve  for  AX,  AY,  AZ  in 
earth-fixed  coordinates.  Also  assume  that  there  w'ere  a  to¬ 
tal  of  m  intervals  during  which  phase  changes  were  measured 
and  the  positions  of  the  satellite  at  the  fiducial  marks  were 
known. 

1  he  solution  of  the  problem  involves  an  iterative 
least  squares  fit  to  t lie  data;  therefore,  initial  estimates  of 
AX,  _>Y,  and  _Z  are  necessary.  From  these,  estimates  of 
-SR  can  also  be  obtained. 

T o  recall  Eq.  ( 12  ): 


SR  =  -  —  (X  -X  +  Y  AY+Z  AZ) 
SR  net  net  net 


Therefore: 

S  (ASR)  _  Xnet 
a  (AX)  "  "  SR 

o  (-SR)  _  ^  net 

S  ( —  Y )  "  '  SR 

5  (wSR)  _  Z net 
5  (-Z)  SR 

For  the  m  intervals,  a  mean  square  error  may  be  defined  by: 

m 

F(_.\, -Y.-Z)  =  L  [x.X  -  (-SR,  -  -SR,  ,)12  . 

k=l  k  k  k-' 

In  order  to  obtain  a  least  squares  fit,  it  is  required  that: 
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3F  3  F  3  F 
3  (AX)  "  3  (AY)  =  3  (AZ) 


To  solve  these  equations,  a  C  matrix  may  be  constructed 
from  the  initial  estimates  of  AX,  AY,  and  AZ;  the  measured 
values  of  N;  and  known  positions  of  the  navigator  and  satel¬ 
lite: 


C(k,  1)  =  NkX  -  (ASRr  -  ASR  ) 
3  (ASR)  3  (ASR) 

r*  (k  9  )  -  _ —  -  ~ 

3  (AX)  3  (AX) 

3  (ASR).  3  (ASR)  , 

ri.  ov  _ k  _ kj^l 

"  3  (AY)  ’  3  (AY) 

3  (ASR)  3  (ASR) 

C(k'  4)  =  3  (AZ)  3  (AZ)  * 


giving  a  C  matrix  of  the  form: 


{C}  = 


C(l,  1) 
C(2,  1) 


C(  1,  2)  C(l,  3)  C(  1,  4) 


C(M,  1)  C(M,  2 )  C(M,  3)  C(M,  4) 


From  the  C  matrix  a  new  matrix  is  defined  whose  terms  are 


A.. 


m 

2 

k  =  l 


9 


where 

i  =  2,  3,  4,  and 
j  =  1.  2,  3,  4. 


This  gives  three  simultaneous  equations: 
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A21  +  A22  d(^X)  +  A23  d(-Y)+  A24  d^Z)  =  0 

A31  +  A32  d<-X)+  As3  d(-Y)+  A^  d(AZ)  =  0 

A41  +  A42  d(-X)+  A43  d("Y)+  A44  d<“Z>  =  0  • 


From  these  three  equations,  d(AX),  d(2Y),  and  d(aZ)  can  be 
found. 


If 

R 

B 

B 

B 

B 

B 

B 

B 

B 


a  new  matrix  is  defined: 


=  A„~ 

A 

A 

11 

33 

44 

13 

34 

21 

=  A_  _ 

A 

A 

A 

23 

44 

43 

24 

=  A 

A 

A 

A 

31 

23 

34 

33 

24 

=  A 

A 

■  A 

A 

12 

32 

44 

42 

!  34 

=  A 

A 

A 

A 

22 

22 

44 

42 

24 

=  A 

A 

A 

A 

32 

22 

34 

32 

2  4 

=  A__ 

A 

A 

A 

13 

32 

43 

42 

33 

=  A- _ 

A 

A 

A 

23 

22 

43 

42 

23 

=  A_„ 

A 

A 

A 

33 

22 

33 

32 

2  3 

and: 


then: 


E=A22,,irA32B21  +  A42n31  , 


dLX)  =  E  (-A21  Bn  +  Agi  B21  -  A41  B31) 


d(”’i  )  =  E  (A21  B12  *  A31  R22  +  A41  B32) 


d(-Z)  =  E  <-A21  B13  +  Agi  B23  -  A41  B33) 
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The  estimates  of  AX,  AY,  and  AZ  are  then  updated: 


AX 

=  AX 

+  d(AX) 

new 

old 

AY 

=  AY 

+  d(AY) 

new 

old 

AZ 

new 

=  AZ 

old 

+  d(A  Z ). 

It  ld(AX)|,  ld(AY)|,  and  ,d(AZ)|  are  all  less  than  a  predeter¬ 
mined  breakout  criterion,  then: 

AX,.  =  AX 
fix  new 

-Y,.  =  AY 

fix  new 

AZ,.  =  ^Z 
fix  new 

Otherwise,  the  new  values  of  AX,  AY,  and  AZ  are  used  as 
improved  estimates  and  the  iteration  process  is  repeated  un¬ 
til  the  breakout  criterion  is  satisfied. 

Now  assume  that  the  antenna  separation  D  is  known 
by  some  other  method  (such  as  surveying)  to  a  high  degree  of 
accuracy.  From  the  relationship: 

D  =  [(AX)2  +  (AY)2  +  (AZ)2]*  , 

it  can  be  seen  that  only  two  independent  variables  exist  in 
this  solution.  As  an  example,  assume  that  AY  and  AZ  will 
be  found.  Then: 


=  [D2  -  (AY)2  -  (AZ)2]* 

ASR  =  -  !  X  [D2-  (AY)2-  (AZ)2]*  +  Y  AY  +  Z  + 

SR  (  net  net  net 
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|OSR)  J_  Xnet  OT) 

3UY)"SR  net‘[D2_  (,y)2.  J^SjJ 


3(1SR)  1  I ,  Xnet  (iZ) 

a<-z)  ”SRI 


The  Taylor  series  expansion  of  FOY,  1Z)  then  becomes: 


1 l 1 

~Y’  «>■£  \  -  <_sV 


iMk-tSRk.i 

d(LY) 


a(:SRk--SRk.1),_j2 
d(_Z)  d(AZ)J  , 

and  the  solution  proceeds  in  a  similar  manner: 


C(k,  1)  =  X.  X  -  ('SR  -  ;SR,  ,) 
k  k- 1 


C(k,  2 )  = 


C(k, 3 )  = 


tJ(LSR)k  d(^SR) 
b(LY)  Ml?) 

a(~SR)k  a(ASR)k_1 

ITZzj  dlZzj 


where 


A  =  E  C  C 
ij  k  =  1  ki  kj 


i  =  Z,  3,  and 
j  "  1.  2,  3. 

A2l  +  A22d<iY,+  A23  d(lZ)=° 


A„.  +  A —  d(AY)+  A..  d(iZ>  =  0 


TMl  JOHNS  HOSKINS  ONIVCSSITT 

APPLIED  PHYSICS  laboratory 

•h.v*»  Spuing  mamtlamo 


d(AY)  = 


d(AZ)  = 


(A31  A23  ~  A21  A33^ 
(A22  A33  '  A32  A23) 

(A21  A32  '  A31  A22^ 
(A22  A33  '  A32  A23J 


AY  =  AY  .  +  d(AY) 
new  old 

AZ  =  AZ  .  ,  +  d(AZ)  . 
new  old 


If  l  d ( A Y )  or  |d(AZ)|  are  greater  than  the  predetermined 
breakout  criterion,  then  the  iteration  is  repeated.  Other¬ 
wise, 


AZ.. 

=  AZ 

fix 

new 

AY,. 

=  AY 

fix 

new 

-X.. 

fix 

=  [D2  -  (AY) 

fix 


AZIMUTH  DETERMINATION 
Recalling  Eq.  (16): 

-SR  =  -  [(X  cos  A  +  Y  sin  A)  cos  <p  +  Z  sin  <pl 
npt  net  1 


net 


-H 

SR 


D 

,  plane  r 

SR  lL(Xnet  C0S  X  +  Ynet  sin  X)  sin  °  ~  znet  cos  oj 


toj  cos  6 


+  [X  sin  A  -  Y  cos  A]  sin  6} 
n**t  net 


From  the  expression  for  -SR,  the  partials  with  re¬ 
spect  to  Dp[ane  and  0  can  be  computed: 
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1  (-'->10  1  ....  .... 

rnr — >  ™  ll<\« cos  x  *  \,et Sl" x)  ° 


piano 


+  IX  sin  X  -  Y  eosX]  sin  0}  » 
net  net  1 


I) 


d  (^SH )  plane  ,  .  .  ,,  .  i  a 

-  77TT  — “rr lL\  sin  X  -  ^  cos  XJ  COS  0 

d  (0)  SR  net  net 

f(X  cos  X  +  Y  sin  X)  sin  ip  -  Z 
net  net  net 


If  we  make  the  following  definition: 


C(k,  1)  =  \.  A  -  (_SR.  -  -SR.  ,) 
k  k  k- 1 


C(k,  2) 


a(-SR)  a(^SR) 


C(k, 3)  = 


o (D  .  )  3(D  .  ) 

plane  plane 

3(_SR)k  9<-SR)kl 


3(0) 


3(0) 


it  follows  that: 


m 


A..  =  L  C  •  C..  , 

13  k  =  1  ki  kj 


where 


i  =  2,  3 
j  =  1.  2,  3 

A21  +  A22d,Vne)+A23  d(9)  =  ° 
AS1+A32  d,Vne)+  A33ti(9)  =  0 


Solving  for  d(D  .  )  and  d(0): 

plane 


Z  cos  <P)  cos  6 
net 

(17) 


cos  <pJ  sin  6}  .  (18) 
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d(D 


plane 


)  = 


^A31  A23  '  A21  A33"' 


^2  A33  "32  "23 


-  A 


d(0)  = 


^A21  A32  '  A31  A22 
^  A22  A33  '  A32  A23^ 


D  .  =  D  ,  .  .  +  d(D  , 

plane  new  plane  old  plane 

0  =  0  ,  +  d  ( 0 )  . 

new  old 


) 


If  ,  d(Dpiane)  or  |d(0)|  are  greater  than  the  predeter¬ 
mined  breakout  criterion,  then  the  iteration  is  repeated. 


Otherwise: 


plane  fix  plane  new 

0,.  =  0 

fix  new 

If  the  antenna  separation  (D)  is  known,  then  the  cal¬ 
culation  is  somewhat  simplified  and  only  the  azimuth  angle  is 
a  variable  of  solution. 

Again  from  Eq.  (18): 

d  (ASR)  Dplane  r  N  N  i  n 

TTeT  SR-  Uxnet  SIn  X  •  Ynet  C0S  COS  9 

-[  (X  cos  X  +  Y  sin  X )  sin  o  -  Z  cos  <pj  sin  0) 
net  net  net 

The  function  F  is  now1  only  dependent  on  0: 
m 

E3  (ASR,  -  ASR  1 )  „ 

[XRX  -  (ASR^  -  ASRR1) - - —  d(0)J  . 

The  C  matrix  then  becomes: 
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C(k,  1)  =  NrX  -  (ASRk  -  ASR  ) 


C(k,2)  = 


3(-SR)  3(:sr) 


a  (6) 


Tie) 


From  which: 


m 


A1  =  k?l  C  k2  ’  Ckl 


m 

A2  =  K?1  Ck22 


A  +  A  d{0)  =0 


Therefore: 


d(6)  =  -A1/A9  . 

The  estimate  of  6  is  then  improved: 

6  =  0  ,  .  +  d(0)  . 

new  old 

If  |  d(0 )  is  less  than  the  breakout  criterion,  then 
®fix  =  0new;  otherwise  the  iteration  repeated. 
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3.  HARDWARE  DEVELOPMENT  AND  INITIAL  TESTS 


INITIAL  TEST  CONFIGURATION 

The  initial  Aztran  tests  employed  two  APL  transloca¬ 
tion  backpacks  as  satellite  receivers.  These  units  (Fig.  11) 
received  the  400  MHz  satellite  signals  on  antennas  separated 
by  about  50  meters.  Both  receivers  were  driven  by  the  same 
5  MHz  local  oscillator  to  provide  phase  coherence  between 
the  two  32  kHz  inr  ^rmediate  frequency  outputs. 

The  phase  difference  between  the  two  32  kHz  outputs 
was  measured  by  the  phase  comparator.  A  flip-flop  was 
turned  on  by  the  negative -going  zero  crossing  of  one  32  kHz 
IF  and  turned  off  by  the  negative-going  zero  crossing  of  the 
other.  The  flip-flop's  output  was  then  integrated  over  many 
cycles  producing  an  analog  voltage  proportional  to  the  phase 
shift  between  the  IF  signals  and  fed  to  a  strip  chart  recorder. 


Results 

A  number  of  satellite  passes  were  received  using  this 
configuration.  The  phase  was  recorded  on  a  strip  chart  and 
the  2  minute  time  marks  were  superimposed  over  the  phase. 
To  determine  the  number  of  phase  counts  during  a  2  minute 
interval,  the  charted  sawtooth  wave  forms  were  counted  be¬ 
tween  the  time  markers  and  recorded. 

From  a  survey,  the  azimuth  of  the  antenna  base  line 
was  determined  and  compared  to  the  computed  results.  The 
errors  in  various  passes  are  plotted  in  Fig.  12  against  the 
number  of  usable  2  minute  data  intervals  received.  As  can 
be  seen,  the  errors  w-ere  strongly  dependent  on  the  number 
of  intervals,  indicating  a  great  deal  of  noise  in  the  data. 

Many  intervals  had  to  be  discarded  because  one  or 
the  other  of  the  receivers  would  occasionally  lose  the  satel¬ 
lite  signal.  If  this  happened  for  even  1  second  then  the  data 
from  the  entire  2  minute  interval  had  to  be  discarded.  This 
initial  test  configuration  demonstrated  the  system  concept' 
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Fig.  1 1  SETUP  FOR  INITIAL  AZTRAN  TESTS 
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however,  it  became  apparent  that  if  tr.o'e  accuracy  were  de¬ 
sired,  better  use  would  have  to  be  made  of  the  available  data. 


IMPROVED  TEST  CONFIGURATION 

There  were  two  basic  problem  areas  in  the  initii  1 
test  configuration.  The  first  problem  was  the  large  amount 
of  usable  data  which  was  being  discarded  when  the  receivers 
would  lose  the  signal  for  a  short  interval.  This  problem  was 
minimized  by  changing  the  measurement  intervals  from  2 
minutes  to  15  seconds.  Thus  only  15  seconds  of  data  would 
be  lost  if  a  receiver  were  to  momentarily  lose  the  signal. 
Shortening  the  measurement  interval  required  a  knowledge 
of  the  satellite  position  every  15  seconds  instead  of  every  2 
minutes.  These  coordinates  were  obtained  by  a  fourth  order 
Lagrangian  interpolation  of  the  2  minute  interval  X-Y-Z  co¬ 
ordinates. 

The  second  problem  involved  data  reduction.  In  the 
original  system,  it  was  necessary  to  manually  read  the  strip 
chart  recording  and  estimate  the  true  phase  from  the  noisy 
phase  recording.  It  became  desirable  to  automate  the  entire 
process,  which  would  allow  data  collection  on  a  24  hour  basis. 
To  solve  this  problem,  use  was  made  of  an  available  Honey¬ 
well  H-21  computer.  This  is  an  8k  computer  with  18-bit 
words  which  was  already  interfaced  to  the  APL  backpack  re¬ 
ceiver  to  perform  the  navigation  calculations. 

A  different  phase  comparator  was  built  that  gave  a 
digital  measure  of  the  phase  difference  in  the  signal.  On 
command  from  the  computer,  one  counter  would  count  the 
number  of  id  MHz  clock  pulses  occurring  in  one  complete 
cycle  of  one  of  the  32  kHz  IF  outputs.  Another  counter  would 
simuuaneously  count  the  number  of  10  MHz  clock  pulses  oc¬ 
curring  between  negative-going  zero  crossings  of  the  two  IF 
outputs.  These  two  counts  were  transferred  to  the  computer 
and  a  ratio  of  them  taken  giving  the  percentage  phase  delay 
to  an  accuracy  of  about  2°. 

The  computer  sampled  the  phase  up  to  100  times  per 
second  and  then  fed  these  samples  through  a  second  order 
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digital  filter.  The  filter  predicted  the  next  phase  measure¬ 
ment  and  the  error  in  the  prediction  was  used  to  generate  a 
correlation  value.  When  a  significant  difference  between  t he 
measured  and  predicted  phase  appeared,  it  was  assumed  that 
one  or  both  of  the  receivers  had  lost  the  signal  and  the  par¬ 
ticular  data  interval  was  disregarded. 

The  value  of  t he  filter's  output  was  stored  every  15 
seconds  and  used  in  the  Aztran  calculations.  Figure  13  is 
a  flow  chart  of  the  steps  involved  in  generating  the  filtered 
phase  output. 

During  a  satellite  pass  the  computer  stores  the  satel¬ 
lite  message  and  the  doppler  counts  used  for  navigation.  Si¬ 
multaneously  the  computer  digitally  filters  the  measured 
phase  and  stores  tne  output  of  the  filter  every  15  seconds. 

At  the  conclusion  of  the  satellite  pass,  the  computer  executes 
the  navigation  calculations  and  prints  the  results.  The  stored 
values  from  the  digital  filter  for  every  15  seconds  during  the 
pass  and  the  satellite's  X-Y-Z  coordinates  are  then  punched 
on  paper  tape  for  later  processing.  When  a  roll  of  paper  tape 
is  accumulated  (usually  in  about  3  days),  it  is  transferred  to 
magnetic  tape  and  used  as  an  input  to  the  Aztran  calculations 
program.  This  PL/I  program  (listed  in  Appendix  A)  computes 
the  azimuth  and  separation,  or  just  the  azimuth  if  t ho  antenna 
separation  is  known. 

Figure  14  shows  Aztran  equipment  installed  in  a  van. 
On  the  left  is  the  teletype  that  prints  the  results  from  the 
computer,  which  is  to  it  right.  In  t he  equipment  rack  (top 
to  bottom)  are  t ho  phase  measuring  interface,  paper  tape 
reader  and  punch,  and  one  of  the  two  satellite  receivers. 

On  the  table  to  the  right  is  a  strip  chart  recorder  and  the 
second  receiver. 

f  igure  15  is  a  photo  of  one  of  the  two  Aztran  antennas. 
It  is  a  10  foot  tall  volute  wound  for  400  MHz.  This  antenna 
was  designed  to  have  uniform  coverage  in  the  upper  hemi¬ 
sphere  .and  attenuate  sharply  any  reflections  from  below  the 
horizon  in  order  to  reduce  phase  errors. 

R«preducMl  from 
b«*l  >viUbU 
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Results 


Tests  were  run  using  the  improved  test  configuration 
for  various  filter  parameters  listed  in  Fig.  13.  From  the 
tests,  it  was  determined  that  the  following  appeared  to  be 
optimum  values: 

K  =  0.  8 

v 

K  =  0.  03  , 

P 

with  a  sample  rate  of  12  samples  per  second. 

With  these  values  fixed,  groups  of  passes  were  re¬ 
ceived  and  analyzed  for  different  values  of  the  constant 
FSHIFT .  Table  1  lists  the  results.  In  part  A,  it  can  be  seen 
that  the  rms  error  in  azimuth  and  separation  is  fairly  con¬ 
stant  for  values  of  FSHIFT  less  than  9.  In  this  calculation, 
azimuth  and  separation  were  both  variables  of  the  solution. 

In  the  azimuth  solution  one  of  the  groups  of  passes  was  ana¬ 
lyzed  solving  only  for  azimuth.  The  separation  was  deter¬ 
mined  by  a  survey.  The  rms  error  in  this  instance  was 
about  half  that  obtained  in  the  azimuth  and  separation  solu¬ 
tion.  This  again  indicates  that  there  is  a  great  deal  of  noise 
in  the  data  since  most  of  these  passes  had  six  or  seven  good 
2  minute  intervals. 
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4.  CONCLUSIONS 


The  tests  run  to  data  have  given  very  encouraging 
initial  results.  The  theory  and  system  concept  have  been 
proven  to  be  practical.  If  the  system  errors  can  be  further 
reduced  to  about  20  of  arc  (a  factor  of  three  improvement), 
then  the  system  can  compete  with  north- seeking  gyro  com- 
passes  in  the  equatorial  and  mid- latitude  regions. 

Multipath  errors  from  ground  reflections  and  track¬ 
ing  filter  errors  appears  to  be  the  dominant  problems.  An 
error  analysis  of  the  Aztran  system  (Appendix  B)  indicates 
that  improved  accuracies  should  be  achievable  if  the  phase 
measurement  errors  can  be  brought  under  control  and  the 
antenna  separation  accurately  measured.  As  the  system 
exists,  without  improvements,  it  is  still  an  unique  tool  for 
polar  navigation  offering  the  only  true  all  weather  azimuth 
measuring  system  which  exhibits  no  singularities  at  or  near 
the  poles. 
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Appendix  A 

AZTRAN  COMPUTER  PROGRAM 


The  following  is  a  listing  of  the  PL/I  program  used 
to  solve  for  antenna  azimuth  and  separation.  The  inputs  to 
the  program  include  navigator's  position,  satellite  position 
as  a  function  of  time,  estimates  of  azimuth  and  antenna  sep¬ 
aration,  and  the  measurements  of  phase  change  during  the 
satellite  pass. 
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PL/I  CPTIHIZ  ING  COMPILES  AZTRAN:  PROCEDURE  OPTIONS  fN  A I K)  RIOREP'; 


SOUICI  LISTING 


STMT  LEV  NT 


1  0  | AZTRAN :  PROCEDUR E  CPTIONS  (MAIN(  RECRDER: 

I/*  AZT  R  A  N  NAVIGATION  PROGRAM*/ 

2  1  0  |  DCL  (ONSOURCE.ONCHARI  UUILTIN; 

3  1  0  |DCL(KNE1(0:69(  ,!NET<0:69|  .  ZN  ET  (0: 6U)  |  FICAT(IS); 

9  1  0  |  DCL  (XS  (C:b4  |  ,IS(D:  EM)  ,ZS  (0:  6<t|  I  FLCAT  (1  5(  s 

5  1  0  (DCL  (SR  <0:69|  ,DSR  (0  :69|  ,DSRD(0  :6U|  ,  DSRAZ  (0:69|  )  t  LOAT(l'): 

6  1  0  |DCL(C(0:63,3’  ,A(2:  3,3(1  PLOAT(15): 

7  1  0  |DCL(IN.IN.ZR,BUTN,BLONG,E,H,G)  FLCAT  (1 5|  S 

B  1  0  |  DCL  PHASE  (0:63| ,  TE£T  (0:63)  ,D  (0:  63)  ,DEL(0:  63|  ; 

9  1  0  I  DC  L  (DA  (7S(  ,TI  (7  S(  ,AZID  (7S|  ,AZIM(75)  )  I  NIT(0)  ; 

10  1  0  (DCL  (SEP(75|  ,INT  (75|  ,  AZIMD(75) ,  AZIRR  (75|  )  I N  I T  (  0)  i 

11  1  0  |  DC  L  NOGCOD  CNAR(80|  .BAD  CHAR  (3(  ; 

12  1  0  |  DC  L  (AZ.AZACT.AZD.AZZR.AZEST.CACT.DAZ)  ELOAT(15); 

13  1  0  |  DCL  (DCPLAN  ,  DELTH.  CEST.DPLAN.P.HT.  R.  RE.  NAVE)  FLOAT  (15)* 

I  /•  •/ 

19  1  0  |  RE-6378190:  /‘EARTH'S  R’DIUS  IN  BETE-S*/ 

15  1  0  |  NAVE*. 79998239;  /*NAV.  "GTh  AT  900BHE.  'BETIBS)  •/ 

16  10  1  P-0. 003353229  ;  /‘EARTH1:.  COEFFICIENT  OF  FLATTENING*/ 

17  1  0  (  BREAK*. 0001;  /‘BREAKOUT  LIBIT*/ 

18  1  0  |  P  = 1 i  /•  PASS  NUMBER*/ 

I  /•  •/ 

19  1  0  ION  ENDPILEITAPE1I  E  EG  I N ; 

20  2  0  I  PUT  SKIP  LISTCEND  CF  FILE'); 

21  201  GO  TO  ND; 

22  2  0  |  END; 

23  1  0  ION  CCNVERSICN  E  EC  I N ; 

29  2  0  (  PUT  SKIP  LIST  (0  NCHAK  .CNSOURCE)  : 

25  2  0  (  DC  1*1  TO  10: 

26  2  1  (  GET  FILE  (TAPE1I  SKIP  EDIT  (NOGOO  D)  (A  (80)  )  | 

27  2  1  I  PUT  SKIP  LIST (NOGCOOI ; 

28  21|  END; 

29  2  0  |  DLA  NK  *0 : 

30  2  0  (  GO  TO  ICAO; 

31  2  0  1  ENO; 

I  /•  •/ 

I  /•  •/ 

32  1  C  .GET  LI  SI (AZACT .A ZEST. DACT . DEST, DELTHI  ; 

I  /*AZACT*SURVETED  AZIMUTH, USE  0  IF  UNKNCNN*/ 

I  /‘AZEST*  ESTIMATED  AZIMUTH*/ 

I  /‘AZACT  t  AZEST  IN  DEGREES  AND  FRACTIONS*/ 

I  /*DACT -SURVEY  ED  SEPARATION,  USE  0  IF  UNKNONN*/ 

I  /‘DEET'ESTIMATED  SEPARATION*/ 

(  /• DELTH* HT  OF  SEflOTE  ANTENNA  ABOVI  RASTER*/ 

(  /‘All  DISTANCES  IN  METERS*/ 

(  /•  •/ 

33  1  0  (GET  lIS1(BLATN,Bt.ONG,HT)  ; 
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PL/1  OPTIMIZING  COMPILE! 


AZTRAN:  PROCEDURE  OPTIONS  (MAIN)  REORDER; 


STMT 

LEV 

IT 

1 

/•BLATN'SURVEIEC  LATITUDE  IN  DECRIES  Or  MASTER  All*/ 

1 

/•BLCNG*SURVEYEC  LONGITUDE  IN  DEGREES  CP  RASTEI  AIT*/ 

1 

/•HT  =  HT  IN  METERS  Of  MASTIB  AIT  ABOVE  GEOID*/ 

34 

1 

0 

1 

GET  LIST  (T H R ESH | ; 

1 

/•  •/ 

35 

1 

0 

1  BE  START :  IS?0 : 

36 

1 

0 

i 

ys*0: 

37 

1 

0 

i 

23  =  0 ; 

36 

1 

0 

i 

ZAr  1; 

39 

1 

0 

i 

BLANK'O; 

i 

/•  •/ 

i 

/•ROUTINE  TO  SIAI'H  fOR  START  OP  PASS  ON  MAO  TAPE*/ 

i 

/•  •/ 

40 

1 

0 

(LOAD*  GIT  MLllTAflll  EDIT  ( E  AC|  (A  1 3|  I  ; 

4  1 

c 

1 

IP  SUPSTR  (BAD.  1,1|  *•  •  THEN  DO; 

42 

1 

1 

1 

6LANK*ELANK*1; 

43 

1 

1 

1 

GET  FILE  (TAPE1I  SNIP  ( 1|  ; 

4u 

1 

1 

1 

GC  TO  LOAD; 

45 

1 

1 

1 

END; 

46 

1 

0 

1 

ELSE  DO; 

47 

1 

1 

1 

If  B I A  A K >7  THEN  DO; 

46 

1 

2 

1 

DDA  T  =  B A D ; 

49 

1 

2 

1 

GET  F ILE  (TAP  E  1|  EOIT  (TTIHE)  (1(31.1  (!))  ; 

50 

1 

2 

1 

BLANK -0; 

51 

1 

2 

1 

GO  TO  OR ; 

52 

1 

2 

1 

END; 

53 

1 

1 

1 

EL S I  DO; 

<•4 

1 

2 

1 

B  LA NK “0 ; 

55 

1 

2 

1 

GET  PILE  (TAPE  1|  SRIF(I)  ; 

56 

1 

2 

1 

GC  TO  LOAD; 

57 

1 

2 

1 

END; 

58 

1 

1 

1 

END  ; 

59 

1 

0 

ION 

DA  (P|  =  D DA  Y ; 

60 

1 

0 

1 

TI  (F| =TTIHE; 

1 

/•  •/ 

1 

/•  READ  IN  THE  IYZS*/ 

1 

/•  •/ 

61 

1 

0 

IDO 

1*0  TO  56  BY  8; 

62 

1 

1 

1 

GET  FILE  ( T  A  P  E  1 1  SKIP  E  C I T  ( I  S  ( I|  ,  Y  S  ( 1 1  ,  IS  ( I|  | 

1 

(E(1  3,  6)  ,1(3)  ,E(13,8(  ,I(3|  ,E(13,8|  |  ; 

63 

1 

1 

1  end; 

1 

/•  •/ 

1 

/•READ  PHASE  All  CORRELATION*/ 

l 

/•  */ 

64 

1 

0 

IDO 

1*0  TO  56; 

65 

1 

1 

|GET 

FILEITAPEII  SKIP  EDIT  (PHASE  (I|  ,D  (I||  (P(7|  , I  (3|  ,P(6;  |  | 

66 

1 

1 

IENC 

; 

1 

/•SCALE  EHASEV 

67 

1 

0 

1 P HAS P»P NASI/ 1000; 

TMt  JOHNS  HOPKINS  UNIVCMSITV 

APPLIED  rnvsics  LABORATORY 

SitVKK  S^M'SG  MaKVLANO 


fi/i  cptibizing  confiler  aztrai:  piocedure  options  <baini  reoboe'; 

SINT  LEV  IT 


68 

1 

0 

1  /‘SEARCH  FOR  LAil  STORED  PHASE  I  ITER  VAL  COD)*/ 

100  1=63  TC  0  81-1: 

69 

1 

1 

1  IF  0(1)  >0  THEN  GO  TO  01; 

70 

1 

1 

1  END ; 

71 

1 

0 

101:  DO  * 8‘TRUNC (1/8 ) ; 

72 

1 

0 

1  /•  •/ 

1  /‘FI  ID  G COD  PHASE  INTERVALS  AID  CCRPUTI  DIFFERENCES*/ 

1  /•  •/ 

1  DO  I«00  TO  2  81  -1 ; 

73 

1 

1 

1  IF  (C|  1|  >THRESH  |D(I)  *0|  THEI  DO; 

74 

1 

2 

1  PHASE  (1|  >999; 

75 

1 

2 

1  ZAP*ZAP* 1; 

76 

1 

2 

1  END ; 

77 

1 

1 

1  ELSE  PHASEdl  -?PASE(I-1| -FHASE(I)  ; 

78 

1 

1 

1  IF  AES  |PHASE(I'|  >15  THEI  PHA  SE  (I)  *999; 

79 

1 

1 

1  ENC; 

80 

1 

0 

1  PH  AS  E | 0 ) =999; 

81 

1 

0 

1 PHASE(1 | =999; 

82 

1 

0 

1 I F  0C-ZAP<2C  THEI  GO  TO  EIOO; 

83 

1 

0 

IIITIPI  «  CC-ZAP ; 

84 

1 

0 

1  /•  •/ 

1  /‘INTERPOLATE  FROH  2  DIN  IITERVALS  TC  15  SEC  IITEiVALSV 

1  /•  •/ 

ICALL  IN1ERPIIS,  00); 

85 

1 

0 

ICALL  IN1EFPITS,  CC)  ; 

86 

1 

0 

ICALL  INTERPIZS,  DO); 

87 

1 

0 

1  /•  •/ 

1  /‘CCFFUTE  NAVIGATOR'S  POSITION*/ 

1  /•  •/ 

1  R*RE/SCRT(COSD(BLATN)“2»(1-F>  “2‘SIID  (BLATI)  “2)  | 

88 

1 

0 

|IN*(R«H1)  ‘COSC  (BLA  Tl|  ‘COSO  (B1CIG)  ; 

89 

1 

0 

III*  (R‘HT)  ‘CCS  0  (BLA1NI  ‘SIID(BLOIG)  ; 

90 

1 

0 

IZI*  |B‘  11- F|  “2»HTI  ‘SI ID  (BLATI)  ; 

91 

1 

0 

1  /•  •/ 

|  /‘CONFUTE  SLAi;  RAIGES  */ 

1  /•  •/ 

IDO  K=0  TC  OC; 

92 

1 

1 

1  INET  1 K |  *IS(R|  -  IN; 

93 

1 

1 

1  T  NET  |  K)  >  TS  (  K|  -IN; 

94 

1 

1 

1  ZIET  IN|  *ZS  (K)  -Zl; 

93 

1 

1 

1  5  R  (R )  ■  SORT  (II ET  (K)  “2»  VIET  (R)  “2‘  ZIET  (R)  “2)  ; 

96 

1 

1 

1  END; 

97 

1 

0 

1  /•  •/ 

1  /‘CONFUTE  THEOI ITICAL  FHASE*/ 

1  /•  •/ 

1  DEL- D; 

98 

1 

0 

|TEST»0; 

99 

1 

0 

|IF  AZACT‘0  THEI  GC  TO  ZZ1 ; 

100 

1 

0 

IIP  CACT‘0  THEI  GO  TO  ZZ1 ; 
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Pl/I  CPTIHIEING  COHPILfR  AZTRAN:  PROCEDURE  OPTIONS  (ntll)  REOICEIi 

STHT  LEI  BT 


101  10 
102  1  0 


10  3  10 

ICu  1  1 

105  1  1 

106  11 

10  7  11 

100  1  1 

1C9  1  0 

110  10 

111  1  1 

112  11 

113  1  1 

11«  11 

115  1  0 

116  10 

111  1  0 

IIS  10 

lio  11 

120  11 

121  1  1 

122  1  1 

123  11 


12U  1  1 
125  1  0 


126  10 
127  1  C 

120  1C 

129  10 

130  10 


131  10 

132  11 


OPLAN'SCRT  <IACT**2-0ELTH**2) ; 

AZ« AZACT ; 

/•  •/ 

/•COrrUTE  PARTI AIS*/ 

/•  •/ 

DC  K  =0  TO  f  C; 

E*-  UNIT  mi  *COSC  (BLONG)  *INET(R)  *SIND  (BLONG)  )  *51  N  D  (  6'  ATI) 
«2N£T(K|  •CCSO(RLATI)  : 

h»  (  (INET  (K|  *CCS  (BLONG)  •»  NET  (R)  *5  I  MD  ( BLOK  G)  I  *COSD(BLATI) 

♦  Z  M  E  T  I  KJ  *SIND|BLATN)  )  i 
G’-INET  (t)  *St  NO  (BLONG)  UNIT  ( K)  •  COS  D  ( BLOI G)  ; 
CSR(R)DPLAN*(f*COSD(AZ)*G*SIND(AZ))/SR(R)  *H*  DEL  TH/SB  (1)1 
END; 

TE  CT  (0(  '991; 

CO  R  =  1  TO  DC; 

T  'ST  IK)  *  |DSR(R-  1|  -DSR  <R) )  /I  Alt  I 
C  EL  |K|  THASE  (K1  -TEST  1R|  ; 

\f  PHASE  IK  |  *999  THEM  DEL(R)-0; 


i zz i :  put  skip  edit ( 1  day  '.oday,1  tire  '.ttire,  'gut  ) 

(A.f  131  .A,E  |U I  .A)  ; 

PUT  SKIf  C>  ; 

PUT  ERIE  WST('  M  AS  IS  -s’  n. PHASE 

T. PHASE  DELTA  PHASE  COIRELATICI  ELEV'); 

DC  R*0  TO  DC; 

Z1»  S u  RT  ( I  N*  *2 ♦  19  •  *2 ♦  ZN* •  2)  ; 

Z2*SvRT(AS(R|**2*1S|R)**2*ZS(R)**2); 

Z3*  (Z2**2-Z1**2-SF(R)**2)/  (2*Z1*SB<R||  ; 

EL=ATANC  (Z3/SQPT  )  1-Z3**2)  3  ; 

PUT  SR  IE  EDIT  (M,IS  IN)  ,TS  <R|  ,  ZS  (I) ,  PHASE  IE),  TEST  (I)  ,DEL,R),D  (I)  ,IL) 
ir  121 .!  119,21  ,P(  19, 21  ,r  (1M, 21  ,rni, *»  .1(31  ,rcB, »),»(*)  ,P(8,«) , 

I  (0|  ,P(8|  ,1(73  ,r(S,1) ) ; 

END; 

PUT  PAGE; 

/•  •/ 

/•  •/ 

/ • SO l V  E  EOF  AZTFUTH  AND  S  E  E  A  RATION*/ 

/•  •/ 

PUT  USTI'AZInuTH  AND  ANTENNA  SEPARATION  ITERATION*  )  ( 

PUT  SAIS  (3|  ; 

AZ* AZEST ; 

DPLAN*SCRT(CEST**2-CELTH**2)  ; 

i«0: 

/•  •/ 

/•COHPUTE  PARTI ALS*/ 

/•  •/ 

START:  DO  R*0  TO  DD; 

E»-  ( IN  IT  (I)  *COSC  (BLONG)  *»NET  (A)  *SIND  (BLONG) )  *SIND  (B.ATI ) 
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tl/l  OPTI II ZI IG  COSPILEI  AZTRAII  PROCEDURE  OPTIOIS  (HAIR)  lEOltf'; 


S TNT  LET  IT 


133 

134 

135 

136 

137 

138 


139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 


152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 


1  1 

1  1 
1  1 
1  1 
1  1 
1  1 


0 

0 

0 


2 

3 

3 

2 

1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 


•  ZNET  (8)  *C  CSD  (BLATN)  ; 

H*  (  (INET  (6)  *COS  D  (B  LONG) *  T  RET  (8)  *SI  ND  ( B LONG) )  ‘COSO  I B’  ATI) 

♦  ZN ET(N) ‘SIND (BLATN) )  i 

G—  I  NET  (6)  *SIN0  (BLONG)  MNET  (8)  •C0SD(BL08G)  ! 

CSBD  (8)  «  I  £*CO  SO  (»Z)  *G*SIN0(AE))  /SR(R)  i 

D5B  A  Z  |  8  *DFLAI*(-E*SIND(AZ)*G*COSD(AZ))/SR(R): 

DS8  IN'  '0  ELAN*  |£  *COS  D  ( A  2)  ♦  G*SIND(AZ))/5B(R)*H*DBLTI/5I(R)  I 
END; 

/•  •/ 

/•CONFUTE  C  8  AT II  I*/ 

/•  •/ 

A*0; 

C*C; 

do  8*0  io  do-1 ; 

IT  PHASI|R«1)  *999  THEN  GO  TO  OUT; 

C  |R,  1)  -  PHASE  (8*1)  *IAVE*DSB  (8*1)  -DSB(R)  F 
0(8,2)  *ISRAZ(8*1)  -CS8AZI8)  ; 

C  (8,3)  *CSRD  1 8  ♦  1 )  -  9S  80  1 8|  ; 

CO  H  *2  TO  3; 

DC  J*1  TO  3; 

A(H,J)*A(H,J)*C|8,H)*C(8,J)  ; 

END  ; 

END; 

OUT:  ENE; 

/•  •/ 

/•CCEPUTE  AND  TEST  RESIDUALS*/ 

/•  •/ 

DAZ*  (A  (2,31  *A  (3, 1|  -  A(2, 1)  *A  (3,3) )  /  (A  (2,  2)*8  13,3) -A  12,3)*A(3,2)  )  t 
D AZ * 5"  '  "*  8*  D  A  Z ;  .  . 

D  D  PI  ,  1)  *A  (2, 2) -A  (2,1 1  *A  (3,2)  )/(A  (2,3)  *»(3,2)  -  A  C,2)*A(3,3) )  | 

DP'  yDOPL  AN ; 

AZ»*. 

If  I>20  THEN  GO  TO  ENOD ;  /»TOC  HART  ITIIITIOIS*/ 

If  ABS  (I  AZ|  >  BRIAR  1HEI  GO  TO  START; 

IF  ABS  IIOPLAN) >0.1  THE*  GC  TC  STIII; 

/*  •/ 

SIGHA-0; 

PUi’eDIK'R  IISIOUAI  ITIIITIOI  RUHB.R' ,  I )  (A  I  r  (  3) )  ! 

DO  8*0  10  DC-1; 

PUT  SR  I P  EDIT(R,C(R,1))  (P(2)  ,I15).P11*»*))! 

SIGHA-SIGHAMC  (8,1)  )  **2; 

If  C  18  •  1)  *0  THE!  HOC  *DDD-1 1 


END; 

SIGHA*S08T  (SIGBA/DDD)  ) 
PUT  SR  IP  ( 2)  ; 

AZD-TRUIC  I  A  2)  ; 

AZIO  (P)  *AZD ; 
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PL/I  CrTIBIZIBG  CCBP IL  £R  AZ'PAN:  PBOCZDURZ  OPTIONS  (BAIN)  BEORCEM; 

SI  BT  LEV  IT 


113  10 

nu  10 

175  10 

176  10 

177  1  0 

178  10 

171  1  0 

1*0  1  0 

i  s  i  io 

1*2  10 

1*3  1  0 


1*7  1  0 

IPS  1C 

1“6  1  1 

1n7  1  J 

1  *  8  12 

1*7  12 

iso  i  i 

171  1  0 

172  1  0 

173  10 

177  1  0 

175  1  C 


176  1  0 

177  1  0 

178  1  0 

177  1  0 

200  1  0 


27  1  10 

202  1  1 

203  1  1 

207  1  1 

205  1  1 

206  1  1 


IAZB=E7*A9S|AZ-AZD|  ; 
lAZIBIi’l  A21; 

ipui  5K : r  eciii'aziputm  ’.kit.’  dig  ',»zb,'  iii', 

I  I  A ,  r  l)|  ,  A,  F  |  7 .  U  I  ,»|  ; 

IPUI  SAIF  EDIT  ('SEPARATION  *,SQBT(CPLAN**2*DELTH**2),a  NE7EBS*) 

I  IA,M7,J|,A|: 

1311  |[|  ';PTIDPLAN**2aDELTH**2l  ; 

ipui  s«:f; 

I ir  A  Z  A  7  T  *  D  THEN  GO  TO  ZZ2  ; 

I  AZSP*60*(AZAOT-AZl  ; 

IPUT  SAIf  EDIT  I  *  A  Z I  RUTH  EBfiOS  *  , A  ZEN  ,  a  BIB')  (A, P  (7, '  )  ,  A)  j 

I  2  Z  2 :  I r  I A  CT  *0  I H  E  A  GO  TO  ZZ7; 

IPUT  5MF  EDITI'SE0.  £  PBO  R  •  , E ACT  -  S  CNT  ( EPL A  I** 2» D1LTH* • 2)  ,  •  BETEBS'I 

I  |A,N7,7|,»|: 

I  /•  •/ 

I  /'TEST  POP  3  S I  DBA  RESIDUALS*/ 

I  /•  •/ 

I ZZ  7 :  J.i  ■  C ; 

ICO  A  7  IC  DC-1; 

IIP  AES  |C  I A .  1)  I  )3*3IGBA  THE!  DO; 

I  PHASE  IA»  1) -977  ; 

I  JJ  JJ»1; 

I  E  N  C  ; 

I  end; 

IPUT  3RIF  ED  IT  |  •  3  Slf.BA  STR  I  P  *  •  ,  3*  S  I  G  B  A )  (A  ,f  (i  ,7 ) )  ; 
l  put  'Alt  121  ; 

IIP  JJ>3  THEN  GO  TO  START; 

iip  c act r o  them  go  to  talli; 

PUT  P  A  U I ; 

/•  •/ 

/•  V 

/•AZINU7H  ITERATION, BUST  A  MOM  DACT*/ 

/*  •/ 

AZTT:  PUT  SNIP  LISTI'AZIBUTH  ITERATICM*)  ; 

1*0; 

PUT  SAIf  |  7|  ; 

DPLAN-3.PTIIACT**2-DLITH**2|  ; 

AZ'AZEST; 

/•  V 

/•COBrilTl  PAPTIAIS*/ 

/•  V 

BEGIN:  10  A  7  TO  DC, 

E  •-  UN  ET  |  A)  *t  CSC  IDLCNGI  «YNET|N)*SIND(BL0NG))*5IID<B-AT») 

•i'Nf  TIA|  *«.'CSt>|BlAlN|  ; 

H  •  |  (INET|A|  *LCS  C(PLOMG)  •  Y  NIT  |A|  *S  I  NO  (ELOMG) )  *COSD(BLATI) 

•ZNET  |  A )  *S  I  ND  I6LATN)  |  *DELTH/CPLAM; 

G*-IMET|A|  •  S  I  NO  |  B  L  0  MG  |  •  Y  N  E  7  |A)*CC50  (  PLOIG)  I 
OSFAZIAI  *DFLAM*  |-E*SIMC(AZ|  •  C*CCS  C  I A  Z  )  ) /S I  (A  ) ; 

DSR  I A )  •DFLAM*  IH*E*COSD  |AZ)  *G*S  I  ID  I A 1 1  )  /S  B  ( A)  ; 


-  OJ 
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T*C  JOHN*  DOMINI  UNIVK*SlTV 

APPLIED  PHYSICS  LABORATORY 

Silvc* 


i 

I 


ft 


Pl/I  CPT I BIZI IG  CC8PIIEK 


illlil:  PIOCEDOIB  0PTI0I5  nun  icoiccr  j 


STHT 

LEV 

IT 

207 

1 

1 

IDS  mm  *tPUI*(H»Z*COSD(IZ)  •  G»  S I  ID  ( A  t>) /SI  (!)  ! 

206 

1 

1 

I  me: 

1  /•  •/ 

|  /»CCHFUTE  C  8 IT 8 1 «•/ 

1  /•  •/ 

209 

1 

0 

|»*0; 

210 

1 

0 

ic*0: 

211 

1 

0 

|  D 0  MO  10  Dt-1  : 

212 

1 

1 

1 1  P  PHIS  I  JK*  1)  *999  THEN  SO  TO  0UTT: 

21 J 

1 

1 

1C  <K.  1)  *FHIS£(K*1)  *6IVE*DSI1  (8*1)  -DSI  (8)  ( 

21» 

1 

1 

|C(K,2)*ISFI7(R*1)-ESFIZ(8); 

215 

1 

1 

1  8  <2, 1)  *1(2, 1)  *C<8,  1|*C  <K,2|  ;  * 

216 

1 

1 

II  (2,21  *1(2,2)  «C(8,2|**2; 

217 

1 

1 

IOU1T:  E ID ; 

1  /•  V 

1  /•COPFUTE  110  TiST  8ESIDUIISV 

1  /•  .  V 

218 

1 

0 

IDI  2  *-57. 29578*1  (2.  1)/l  (2.2)  i 

219 

1 

0 

1 IZ*IZ<  EIZ ; 

220 

1 

0 

1 1  - 1* 1 s 

221 

1 

0 

IIP  I>20  THE!  GO  TO  T!UI; 

222 

1 

0 

IIP  IPS  (C8Z)  7B8EIK  1HEN  GO  TO  tEGIIj 

1  /•  •/ 

223 

1 

0 

|SIG6»*0; 

229 

1 

0 

1  ODD*  DO; 

225 

1 

0 

IPUT  SKIP(I)  ; 

226 

1 

0 

1  PUT  EDIT  ('  R  IZSIDUIL  ITEIITIOI  lUFLEI'  , 

227 

1 

0 

IDO  K  •  0  TO  DD-li 

228 

1 

1 

IPUT  SKIP  EDIT  (K,C(K,1)  )  (M2),l(5),r(19,9))i 

229 

1 

1 

|  S I CHI  *  SIG8I * (C (8 , 1) ) **2  3 

230 

1 

1 

1  IP  C  (K  .1)  *0  THEI  DDD*DDC*1  i 

731 

1 

1 

1  EN  C  i 

232 

1 

0 

ISIGHI'SCM  (SIGHI/DCD)  ; 

233 

1 

0 

IPUT  SKI  F  12)  ; 

23« 

1 

0 

1  IZ  C*TR  UNC  (1  Z|  ; 

235 

1 

0 

1 IZIPID (FI  *IZC; 

236 

1 

c 

|»ZH-60*IbS(IZ-IZD)  ; 

237 

1 

0 

IIZI6HIPI  *IZH; 

238 

1 

c 

IIZEI*60*(IZICT-IZ,  : 

239 

1 

0 

IPUT  SKIP  EDIT  ('IZIHUTH  MU,'  DIG  ',1X8,'  HI’) 

|  (1 ,  P  (3|  .  1 ,  P  (7,9)  ,1)  ; 

290 

1 

0 

IIP  IZICT'O  THEI  GO  TO  ZZ8 ; 

291 

1 

0 

1  /•  •/ 

|  /‘TEST  POI  3  SIG8I  IISIDOIISV 

1  /•  V 

292 

1 

0 

IZZI:  JJ *0; 

293 

1 

c 

|D0  K*0  TO  DD-1; 

299 

1 

1 

IIP  IIS(C(8,  1)1  >3*SIG8I  THEI  DO; 
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T*i  JOHNS  HC"«  ns 

Fr  '  v  ..ABC«  at  'CV 

^  f«  5p*  «.i  -.•fits*' 


Pl/I  CPTXHIZING 

COflPILFR  AZ 

TIER:  PROCEDURE  OPTIONS  { HA  I  R| 

REORDER; 

STRT 

LEV 

IT 

*45 

1 

2 

1  P  H  A  5  f  (M  1)  >999 

2  46 

1 

2 

1  J  J  *  J  .1  ♦  1 ; 

24  7 

1 

I 

1  END: 

244 

1 

1 

J  F  *»  C  ; 

244 

1 

C 

1  P  in  SKIP  t  OXT  ( 1  3 

DIGHA  STPIf  3*SIGHA|  (A, Tib, 4)1  ; 

250 

1 

0 

1  PUT  S* 1 1  Cl  ; 

*51 

1 

0 

1  IF  JJ>C  THE  N  i‘.C  T 

C  EIGIN; 

252 

1 

c 

1 T  ALL  T ;  PP  ♦  1 ; 

25  J 

1 

0 

1  I F  F>75  THEN  P*75 

254 

1 

0 

IENCC:  PUT  PAGE; 

255 

1 

0 

|GO  TO  RESTAFT; 

1  /• 

•/ 

1  /• 

•/ 

1  /*LAGRANGIAN  INTERPOLATION  ROUTINE*/ 

1  /•THE  ORDER  Or 

INT  EPPOLAT ICN  VARIES  WITH  THE 

RUHBEE  OP 

INTERVALS*/ 

1  /* 

•/ 

.''4 

1 

c 

IINTERP!  PROCEDURE  (ARMY,  CC)  PECBDER; 

257 

2 

c 
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Appendix  B 

AZTRAN  ERROR  MODEL 


SUMMARY 

A  mathematical  error  model  of  the  Aztran  system  has 
been  developed  to  study  the  limiting  accuracies  of  the  system. 
Reasonable  measurement  uncertainties  were  selected,  and 
used  as  inputs  to  the  model.  Based  on  these  uncertainties, 
a  50  meter  baseline  system  should  be  able  to  achieve  an  ac¬ 
curacy  of  about  50"  to  60"  of  arc  rms  and  a  100  meter  base¬ 
line  system  should  approach  25"  of  arc  rms  in  the  mid  lati¬ 
tudes.  The  values  vary  somewhat  with  the  azimuth  angle  of 
the  antenna  and  the  pass  geometry. 

Uncertainties  of  20  meters  rms  in  latitude  and  longi¬ 
tude  and  4  meters  in  height  were  used  for  the  navigator's  po¬ 
sition.  This  degree  of  accuracy  is  obtainable  by  averaging 
a  small  group  of  refraction  corrected  passes.  The  error  in 
the  navigator's  position  contributes  almost  nothing  to  the 
total  system  error  and  could  be  10  to  100  times  greater  in 
magnitude  before  it  approached  the  value  of  some  of  the  other 
errors. 


An  rms  value  of  8  meters  in  positional  uncertainty 
was  assumed  for  the  satellite  in  each  of  the  three  Cartesian 
coordinates  at  any  one  fiducial  mark.  This  value  of  8  meters 
is  combined  satellite  orbit  prediction  error  and  H-21  com¬ 
puter  conversion  error  from  Keplerian  to  Cartesian  ele¬ 
ments.  No  error  correlation  was  assumed  from  measure¬ 
ment  to  measurement.  The  position  errors  contribute  little 
to  the  system  error  and  could  be  as  much  as  a  factor  of  10 
larger  before  they  u'ould  begin  to  become  significant. 

Two  different  pairs  of  uncertainties  were  used  for 
estimates  of  the  antenna  separation  and  height  difference. 

For  the  first  pair  3  cm  and  1  cm  were  used,  respectively. 

At  this  measurement  accuracy,  these  uncertainties  tended 
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to  be  the  dominant  error  source  for  a  100  meter  baseline 
system. 

The  second  set  of  uncertainties  were  1  cm  in  separa¬ 
tion  and  0.  5  cm  in  height  difference.  These  values  require 
considerably  more  precise  measurement  techniques  but  are 
still  feasible.  In  this  case  the  antenna  measurement  errors 
were  significant  for  a  50  meter  baseline  system,  however  in 
a  100  meter  baseline  system  the  dominant  error  source  was 
the  phase  measurement  error. 

In  the  current  Aztran  system  the  phase  is  sampled 
100  times  per  second  by  the  computer  in  a  second  order  dig¬ 
ital  filter.  The  maximum  phase  frequency  being  followed  is 
about  1  Hz.  Sources  of  phase  error  are  jitter  in  the  receiver 
phase-locked  loops,  multipath  reflections,  and  sampling  error. 
A  value  of  0.01  cycle  of  phase  error  was  selected  as  the  rms 
measurement  error.  This  is  not  to  imply  an  error  of  0,01 
cycle  per  phase  sample,  but  that  the  smoothed  output  from 
the  digital  filter  (which  is  stored  every  15  seconds)  will  have 
this  rms  error.  Due  to  the  high  sampling  rate  (1500  samples 
between  stored  values),  the  error  from  stored  value  to  stored 
value  was  assumed  to  be  uncorrelated. 

The  second  section  is  a  development  of  the  equations 
used  in  the  error  model  and  the  final  section  discusses  the 
results  of  the  model  in  more  detail  and  their  implications. 


ERROR  ANALYSIS  EQUATIONS 

In  the  Aztran  system,  a  signal  path  length  difference 
exists  between  the  two  antennas.  This  difference  is  deter¬ 
mined  by  measuring  the  phase  difference  in  the  two  received 
signals: 

-SR  =  NX,  (B-l) 

where: 
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ASR  is  the  slant  range  difference, 

N  is  the  number  of  wavelengths,  and 
^  is  the  wavelength  of  the  received  signal. 

Recalling  Eq.  (12)  of  this  report: 

ASR  =  (X-AX  +  Y-AY+  Z-AZ), 

where,  in  Cartesian  coordinates: 

X  =  X  +  ...4  -X  .  =  X  -  X 

satellite  navigator  s  n 

Y  =  Y  -Y  =  y  -  Y 

satellite  navigator  s  n 

z  =  z  t  -  z  .  =  z  -  z  . 

satellite  navigator  s  n 
From  Eq.  (10)  of  this  report: 

D  =  (AX2  +  AY2  +  AZ2)^  , 


where: 


D  =  the  antenna  separation  and 
^plane 


=-J D2  -  (AH)2 


(B-2) 


(B-3) 


From  Eqs.  (13)  through  (15)  of  this  report: 

AX  =  -  D  .  [ sin  <p  cos  X  cos  0  +  sin  X  sin  0] 

plane 


+  _  H  c  os  <p  c  os  X 


(B-4) 


:Y  =  -  D  .  [ain<p  sin  X  cos  0  -  cos  X  sin  0] 
plane 


+  AH  cos  <p  sin  X 


,Z  =  D  ,  cos  <p  cos  0  +  AH  sin  o  , 
plane 
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v here: 


<P  =  navigator's  latitude, 

X  =  navigator's  longitude, 

0  =  baseline  azimuth  angle,  and 
-H  =  antenna  height  difference. 

From  the  preceding  equations,  iSR  can  be  defined  as: 

-SR  =  f(o,  X,  HT,  X  ,  Y  ,  Z  ,  D  ,  ,  iH,  0)  (B-7) 

s  s  s  plane 

The  true  J>SR  can  be  expressed  in  terms  of  a  Taylor  expan¬ 
sion  about  the  nominal  value: 


-SR  =  _SR  +  ~  (O  -  <P  )  +||  (X  -  X  )  + 
t  n  t  n  ^X  t  n 


=  £SR  +  |--6o +  |y6X  +  ^y-6HT  + 
n  do  dX  OUT 


The  measured  ..SR  can  be  expressed  in  terms  of  the  true 
value  and  the  measured  phase  error: 


iSR  =  ASR  +  —  6  (ph  ) 
m  n  360  m 


By  subtracting  the  nominal  value  from  the  measured  value, 
the  true  value  drops  out: 


SK 


S|< 


A 

3  00 


6(p!i  ) 
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a  e 
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During  a  satellite  pass,  many  measurements  are 
made.  An  estimated  azimuth  angle  is  then  improved  in  an 
iterative  least  squares  sense  to  drive  the  nominal  £SR  to  the 
same  value  as  the  measured  A SR.  Equation  (B  -8)  now  be¬ 
comes: 


af 

ae 


60 


^rSiph  >-[|£ 

360  m  o<p 


af  . 

aAH_ 


6<p 

6X 


(B-9) 


|_6AhJ 

In  a  similar  manner,  the  antenna  separation  error  sensitiv 
ity  can  be  determined: 


I 

E 

I 

I 
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—  6D  = 

aD  360 


t/  ,  \  r  3  f  3f 

6(phm’-LSsT 


af  af.. 
'  ’  aaH  ae J 


69 

6X 


03-10) 


6^H 

_60  _ 

Assume  that  during  the  satellite  pass,  n  separate 
measurements  are  made,  each  representing  a  data  interval. 
Eq.  (B-9)  now  takes  the  form 


Sfl 

»f, 

9f. 

S(phm>l 
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1  ,  , 

ax 

1 

8j|( 

Sf2 

8  0 

60  =  360 

<1(phm,2 

- 

5f2 
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5f2. 

8X 
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Let  Eq.  (B-ll)  be  represented  in  matrix  notation  as: 

—  *  60  =  3?0  *  —  '  —  '  —  *  (B-12) 

Using  Eq.  (B-12)  and  estimates  of  the  various  errors  ( O ,  X, 

.  .  .  .,  -iH,  6ph  )  an  error  in  60  can  be  determined  to  be 
m 

69  =£  '  (no'  “'S'  er)  •  (4-J  •  <B'13> 

Houev’er,  Ahe  desired  result  is  not  the  error  from  a 
single  pass  for  a  given  set  of  biases  but  the  statistical  vari¬ 
ance  for  a  pass  given  an  estimated  variance  in  the  error 
sources. 


This  becomes: 


•5<h'  -  K! 


+  K  ' 


1 1: h  ■  rn‘ . 


r  ■ 


(B-14) 


In  Eq.  (B-14)  [M_'  ]  is  the  measurement  error 

variance  matrix  which  is  diagonalized,  implying  no  correla¬ 
tion  from  measurement  to  measurement.  Similarly  [ER  • 

ER:  ]  is  the  diagonalized  error  variance  measurement,  again 
assuming  no  correlation  from  error  source  to  error  source. 

From  Eq.  (B-14)  the  system  sensitivities  to  each  var¬ 
iable  can  be  derived  and  a  total  error  budget  can  be  deter¬ 
mined.  The  only  remaining  problem  is  to  calculate  the  par¬ 
tial  derivatives  of  f.  Recalling  Eq.  (B-2): 


-SR  =  -  —  •  (X  •  IX  +  Y  •  1Y  +  Z  •  iZ) 
oK 


Therefore,  the  partial  derivative  of  1SR  can  be  written  as: 


3  SH 


Ml  /SR  > 
5  i 


<  \  -\  +  \  IV  +  7.  17.) 


y  .  y  +  /.  +y.i=Z)  .  (B-  15) 

dj  3i  3 i  31  I 
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where  i  is  one  of  the  variables  O,  A  ,  Xg,  Yg,  Zg, 
Recalling  Eq.  (9)  from  this  report: 


1/SR  =  (X2  +  Y2  +  Z2)'* 


so  that: 

3  (1/SR) 
3  i 


srl 


S*v£  +  zf 


r) 


X  =  X  -  X 
s  n 


J 


X  =  X  -  (R  +  H  )  cos  O  cos  A  , 
s  1 


where: 

R  =  Re  •  [cos2  O  +  (1-f)2  sin2  oJ  ^  , 

Re  =  equatorial  earth  radius, 

O  =  navigator's  latitude, 

A  =  navigator's  longitude,  and 

H,j,  =  navigator's  height  above  the  reference 

From  Eq.  (B-17)  it  can  be  seen  that: 


s 


■r~—  =  -  COS  O  COS  A 

0  H^. 

3X 

g^-  =  (R+HT)  sin  0  cos  A 

R3  2 

+  — =  •  f  •  (f- 2 )  •  cos  O  sin  O  cos  A 

Re 


,  ah,  e. 


,  (B- 16) 


(B-17) 


ellips  oid. 
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=  (R  +  H  )  coso  sin  X  . 
0  A  1 


V  =  Y  -  Y 
s  n 


Y  =  Y  -  (R  +  H  )  coso  sin  X 
s  1 


(B-  18) 


Therefore: 
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d  ^  1 

3  2 

+  R  •  f  *  (f-2)  •  cos  o  sino  sin  X 
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a_v 

ax 


(R  +  H  )  cos  o  cos  X  . 


Z  =  Z  -  Z 
s  n 


Z  =  Z  -  [R(l-f)2  +  H„]  sinO 

s  T 


(B- 19) 
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!=•  =  -  [R(l-f)  +  H_ ]  cos  to 

0(0  1 


R  2  2 

+  .  (1-f)  .  f.  (f-2)  sin  (0  cos  O 

Re2 


AX  =  -  D  (sin  (3  cos  A  cos  0  +  sin  A  sin  6) 

plane 


+  AH  cos  <0  cos  A 


(B-20 


=  -  sin 


to  cos  A  cos  0  -  sin  A  sin  0 


plane 


a_X  D  (sin  (0  cos  X  sin  0  -  sin  X  cos  0) 

ft*  pla"e 


cos  to  cos  A. 


d ^  =  -  D  cos  >0  cos  X  cos  0  -  -H  sin  O  cos  X 

do  plane 


=  p  (sin  0  sin  A  cos  0  -  cos  A  sin  0) 

9  A  plane 


-  _H  cos  o  sin  A  . 


’  Y  =  -  D  (  sin  O  sin  A  cos  0  -  cos  A  sin  0) 
plane 


+  AH  cos  p  sin  A  . 


(B-2 


-£±L_  =  -  sin  o  sin  A  cos  0  +  cos  A  sin  0 


plane 
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plane 

-  ■jJl  sin  P  sin  X 


-13  (sin  P  ccs  X  cos  0  +  sin  X  sin  0) 
plane 

+  MI  cos  P  cos  X  . 
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plane 
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RESULTS 


The  results  of  th;  error  analysis  for  various  passes 
has  been  tabulated  to  show  the  total  rms  error  and  the  con¬ 
tribution  from  each  error  source  (in  seconds  of  arc)  for 
antenna  arimuths  from  0°  to  90°.  The  pass  geometry  is  de¬ 
scribed  for  each  pass.  The  error  values  chosen  were  2  0 
meters  in  latitude  and  longitude,  4  meters  in  height,  8  meters 
in  the  satellite's  X,  Y,  and  Z  coordinates  and  0.01  X  in  mea¬ 
surement  noise. 


Table  B- 1  shows  the  results  of  a  few  passes  over 
APL's  Howard  County  facility  with  a  100  meter  baseline  for 
antenna  measurement  errors  of  1  cm  in  separation  and  0.  5 
cm  in  height  difference.  Table  B-2  shows  the  same  passes 
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over  APL  with  a  50  meter  baseline  and  antenna  measurement 
errors  of  1  cm  in  separation  and  0.5  cm  in  height  difference. 
As  can  be  seen  by  comparing  Tables  B- 1  and  B-2,  azimuth 
error  due  to  errors  in  antenna  separation,  height  difference, 
and  phase  measurement  vary  inversely  with  antenna  separa¬ 
tion,  and  the  others  remain  constant. 

Figures  B- 1  through  B-4  show  the  rms  azimuth  error 
in  seconds  for  various  antenna  azimuths  and  pass  elevations. 
From  these  graphs  it  can  be  seen  that  the  key  to  system 
improvement  lies  in  decreasing  the  antenna  separation  mea¬ 
surement  error  to  0.  1  cm  or  less.  This  accuracy  requires 
laser  ranging  but  is  certainly  possible. 

A  listing  of  the  computer  simulation  used  to  generate 
the  data  in  Tables  B-l  and  B-2  is  given  in  Fig.  B-5. 


TABLE  B  1 

RESULTS  OF  PASSES  WITH  100  METER  BASELINE 
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TOTAL  ERROR 


—  SEPARATION  ERROR  (1  cml 
—  —  PHASE  ERROR  (0  01  X) 
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2  i- 


1 _ J _  1 _ I - 1 - 1 - 

0  30  60  90  60  30  0 

E  —  |  »-  W 

SATELLITE  ELEVATION  AT  CLOSEST  APPROACH  (degrees) 

Fiq  B  1  RMS  AZIMUTH  ERROR  FOR  SATELLITE  PASSES  OVER  APL  (AZ  =  0°) 
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Fig  B  2  RMS  AZIMUTH  ERROR  FOR  SATELLITE  PASSES  OVER  APL  (AZ  =  30°) 
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HMS  AZIMUTH  ERROR  Imc  seconds) 
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Ft)  B  3  RVIS  AZIMUTH  ERROR  FOR  SATELLITE  PASSES  OVER  APL  |AZ  =  60  ) 


RMS  AZIMUTH  ERROR  (arc  seconds) 
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Fig.  B-4  RMS  AZIMUTH  ERROR  FOR  SATELLITE  PASSES  OVER  APL  (AZ  =  90°) 


Pl/I  OPT  I A IZ I NG  CORPILtS 


AZ  A  L:  PROCEDURE  OPTIOR5  (AAIN)l 


SOURCE  LISTING 

STRT  LET  NT 


1  0  IAIAL:  PROCEDURE  CPTIORS  (BAIN)  I 


2  1  0 

3  1  0 

u  1  0 

'  1  0 

6  1  D 

7  1  0 

A  1  0 

9  1  0 

10  1  0 

11  10 

12  10 

13  10 

1«  10 

15  1  0 

16  10 


17  10 

1  A  10 

19  10 

20  1  0 

21  1C 

22  1  0 

23  10 

24  1  C 

25  1  0 

26  1  0 

27  1  0 

2A  1  0 

29  1  0 

30  1  0 


31  1C 

32  1  0 

3  3  10 

14  1  0 

3 e  1  0 

36  ID 

37  1  0 

3«  1  0 

39  1  D 


DOL  IS  (A)  ,  TS(8|  ,ZS  (0)  ; 

DC  L  PI(A|  ,PKA|  ,PZ(8)  i 
DC  L  rDI  (A|  ,  PD  1  (B|  ,PDZ  (8)  ; 

DCL  FR  ( A |  ,5AZ  <A|  ; 

DC  L  PR  <A|  j 

DCL  MS,  A|  PL  CAT  ; 

DC l  PER  <8,A|  FLOAT; 

PCI  RTF  (8|  ; 

DCL  RTPER  (8)  ; 

DCL  LAT  FLOAT ; 

DCL  LOR  FLOAT  ; 

DCL  A  FLOAT; 

DOL  AT?  FLOAT; 

DCL  N  (2|  CHAR  (81  INIT  ('NCNTH*  ,  'SOUTH*)  ; 

DCL  E  (  2|  CHAR  (4)  INIT(*EAST',  •  NEST  •)  ; 

/•  CONSTANTS  •/ 

PE  =6)78140; 

RAVE--.  74948234;  /*NAf ELEHGTH  AT  400HHZ.  •/ 

F  *  1/2 98 . 22  ; 

PP  =  (1-F>**2; 

DTP-ATAN  (1J  /AT  AND  (1)  ;  /*DEG  TO  RADIAIIS  •/ 

AT  R  =  DT A/60 ;  /•  AIR  TC  RADIANS  •/ 

RTS-60/ATR;  /•  RADIANS  TO  SECONDS*/ 

DCL  EL  (0 : 9)  INIT<-1, 0,0, 0,0, 0,0,0. 0,-1)  ; 

pi-0: 
pi  =0; 

PZ-0; 

PDI’O; 

PD  1*0; 

PCZ-0; 

/•  INPUT  DATA  •/ 

GET  II  AT (LAT, ICN.HT)  ; 

GET  SP IP  LIST  <D,DH|  ; 

GET  SNIP  LIST (£R,A| ; 

PUT  TP  IP  EDIT!'  ','LAT','LOR','HT','IS  •  ,  •  IS  •  , 'ZS  '  ,  •  D' ,  •  OH*  ,  •  ■• ) 
(I  ( 111  ,A  ,  I  (6|  ,2  (A  ,1  (A)  |  ,*  (A,  I  (9|)  ,A,  1(10)  ,A,I(9)  ,A)  ; 

PUT  SKIP  F  C  IT  ( ER  ,  A|  (I  ( 19|  ,  2(E  (7,  3)  ,  K  4)  )  ,6(P(5,3)  ,1(6)  )  ,F  (4,2))  ; 

put  sp  IP  (4)  ; 

PUT  ED:T('LAT*',TRUNC(LAT|  , 'CEG.'.OO^ABSITRONC  (LAT) -LAT)) 

(I  ( 10|  ,  A  ,P  (4|  ,1  (1)  ,A,F  (8,4)  ,  I  (1)  ,  A)  ; 
put  SPIPdl  ; 

PUT  EDIT  ('ION*',TNUNC<LON)  , • DBG. • , 60* ABS (TNONC  (LON)-IOB)) 


F  g  B  5  LISTING  OF  AZTRAN  ERROR  ANALYSIS  SIMULATION 
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THl  JOHNS  KO»»  NS  UNIVERSITY 

APPLIED  PhyS'CS  LABORATORY 


S'LYfH  S#«'N<1  Mi 


ft /i  oi>  rim  line  coKfiim 


UiL:  PPOCIDOPI  CPTIOPS  (Mill  I 


3T«T  LIT  »T 


1  1  I 
1  1  I 
1  1  I 


II  110)  ,i,r  |«)  ,1  (1|  ,i.f  (8,U|  .1(11  ,i|  ; 

PUT  SUP  tOIT|'HIICHT.*.HT.*«. '1(1(101. *.K7.2),I(1|,*J, 
PUT  SUP  IDIT  (*D*‘  ,0,  •  II.  •  |  (I  |10)  ,l,P(t.2)  ,i|  ; 

PUT  SUP  I0IT  I # D H* *  .  DP.  •  *|  ( I  1 10|  ,1,  P  («,  ",  ,1)  • 

PUT  PAGE;  ' 

tit  in  »oi  i)  /ms 

I*  12)  ■  IP  |2)  /  (COS  D  (LIT)  *  PE)  | 

P*P*PATE/)60i 

/•  */ 

I IT: DO  1*1  TO  «j 

GIT  SUP  LIST  (I3|I)  .TS  |I|  .13(1))  | 

EM: 

CLl*COSD  IL1T)  ; 

CLC*COSD  |  LON)  ; 

S  Li'*  SI  ND  (LIT)  i 
SLO*SIPD|LOP)  ; 

p-pe/sqpt  icLi**2«rr*sii**2| : 

II  -  (P*  NT)  *CL1  *CLO ; 
m«l*PHT)  •CIMSIO: 

Z9* (P*II *HT) • SLi; 

PI  1 1)  ■  (S  »HT|  *SL1  *CLO*P**3*P*  (  P-2)  *Cll**2*SLA*CLO/»E**2  | 
PI (2) • TP; 

PI  13)— Cll*CLO: 

PI(«)  «1j 

PI  (1 )  •  |P  ♦HT|  • SLl*SLO«l**]*P*  (P-3)  •CLI**2*Sll*SLO/PI**2  « 
PI  12)  ■-  IP  ; 

PI  |3)«-Cll«SLC| 

PT(S)*1: 

p:  ii)«<>(P«rr*iTi  •cLi«i**]*rr«r*(r-2)  *sla**2*cla/»i**2: 

PI  (  3)  ■•SLI ; 

P  2  (6)  *1 J 

0O  IZ-0  TO  90  IT  10 ; 
sii*0: 

P  *0; 
pp»0: 

Cl  *COSD ( I Z)  ; 

S1»SIPD(1Z)  ; 

DI*-D*  (SLl*CLO*Cl ♦  SLO*Sl)*DH*CLI»ClOi 
DT  —  D»  |Sll*SLC»Cl-ClO*Sl)  »DI*CLI*3LO; 

DZ»D*CII*CI*DP*SLli 

/•  COPPUTE  PUT I1LS  •/ 
PDIl*D*(Sll*CLO*Sl-Sll*CI|  ; 

PDT1*D*  (3ll*SlC*Sl*ClO*CA)  t 
PDZI—  D*CLI*SI; 

PDIdl  —  D*(Cll*ClO*Cl)  *0*311  *CLO| 

PDI  (2)  *-01: 
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’  -t  .»>•■** 
c? 


•  -  ABT  -  -  ”  ' 


'  r* 


fisi  cpTimnir  ccmhleb 


azai:  ptoccooai  omois  (maip)  i 


st«:  Lff  i: 


•  )  i  *, 

**4  i  i 

•5  1  1 

P6  1  1 

«7  1  • 

«•»  1  1 

“6  1  1 

70  1  , 

61  11 

1  1 

61  12 

6*  i  ; 

65  1  2 

66  12 


P0I*7»  - «IA*C IC*CA-SLO*SA; 

PDI  (0) ■CL»*CtC; 

POIfl)  «-t*C LA*SLO»CA-DN*SLA*SLO; 

PDI  (2) *  Df ; 

PD!|7»  *- eLA*SlC*CA*CLO*3A; 

PDf  (*•!  *CIA*SLC; 

PDZ(I)  *- D*  SLA •CA*DH*CLA ; 

PD7(7>  CLA*CA  ; 

PCZ  |0)  * SLA; 

DO  1*1  TO  6; 

I  -  IS  f 1 1  -IN; 
f  MSI  I)  -»«; 

Z*Z5(!) -ZP; 

sa*S;PT  T*»2*Z**2)  ; 


77  1  J 

6“  1  2 

<*<»  1  2 

i:o  i  z 


/•  COMPOTE  ELEV ATIOP  •/ 
A  4 *1 N*T* Y 4* Z*Z»; 

S*  rp*pr*«2/p  ; 

SLID  ATAur  'A *9/ (B**2-A**2))  l 
IP  rl(  ;»  <0  TMf  P  r.O  TO  A  A  ; 


m  ir 

*C2  i  2 

10)  1  2 

104  1  ; 

US  1  ; 

1C6  1  ; 

107  1  2 

i"6  1  2 

106  i  j 

116  1  ) 

111  1  ) 

112  1  * 

H>  1? 

114  1  2 


/•  coarurr  azimuth  •/ 
a  »{•  if  * « r  i  •  m-is  ti)  mp> ; 

*»  *5 HI  •PE«6  2-ZW*(«S|I)4III#TS|I)»f»*XS|I)  MU  j 
SAZC)  ATA4CCA/P)  ; 

A  7  (!)  /SjPT  IIS  CI)««2*TS(I)  ••2*XS  (X)  •  •E)  ; 

3  A05  (SLA)  ; 

T  P  A  <  i*  --TB  'AZ(t)  >5At(X)#1fO. 

IP  jAIIIkO  T98A  SAZ  (X)  >SAZ(X)  •  360  ; 

■P  C) T -( I»PD« A*I*PDT A*1*PD? )  /S  P; 

30  J»1  TO  6; 

A  (l*3l*T*DT*Z*PZ)  • ( I • P I «J)  *T*PV(J) »  Z*  P Z  (J)  ) /S  •••  )  ; 

*  •°,#T*fCf  <J»  *rr  W)6Df*f6tOZ  (J)  *P*(J)  •0I)/SP; 

P6D;  /•  ftP  J  •/ 

AA:  l«0;  /•  EMC  I  •/ 


115  1 

11*  1 
117  1 

»14  1 

116  1 
120  1 
121  1 
122  1 


/•COMPUTE  PASS  LIMITS  (CAN)  •/ 

30  I  4  TO  0  AT  - 1 ; 

IP  U(X)  <0  THPt  00  TO 

f4C; 

PB:  0*1*1; 

30  I  S  TO  6 ; 

IP  81(1)  <0  TMf  ■  CO  TO  CC; 

BAD; 

CC:  CM; 


12) 


PTP-O; 


/•COMPOTE  ft  TMAaSPOSE  MV 
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TNI  JOhNI  KO*»  N*  UNIVt*5ITV 

4PFL  Er  P'-VS'CS  LABORAT CRV 

V-lv“  w«au*xo 


tl/1  CPTIHI1IEG  COMPILES 


ItlLl  EEOC  EDO  El  CPTI0E3  (SHE)  t 


STMT  LEV  ST 


12V 

1 

1  i 

00  1*0  TO  H; 

125 

1 

2  1 

ET8'ETE‘EE  ( I)  **2i 

126 

1 

2  | 

EMC; 

1 

1 

/•  COMPUTE  ST*  •/ 

121 

1 

1  | 

8TX “0; 

128 

1 

1  | 

DO  J-1  TO  S  ; 

129 

1 

2  | 

00  I-G  TO  Hi 

1  V 

1 

3  1 

9TMJ)  *ETX(J)  •  El  ( I)  •«(.",  J)  i 

131 

1 

3  1 

ESC; 

13k 

1 

2  | 

EH  0  * 

1 

1 

/•  COMPUTE  I  EE  EET  XT  •/ 

13  3 

1 

I  | 

XEE-0; 

13V 

1 

1  | 

oo  j-g  to  hj 

115 

1 

2  1 

00  X  -G  TO  Hi 

1  16 

1 

3  1 

00  IH0E1-1  TO  8; 

111 

1 

V  | 

KEH(J.I)  •XES(J.I)  ‘X  ( J#  IS  DEI)  vEB|XE0EX)**2*X|IvIEDEE){ 

119 

1 

v  1 

EMC: 

119 

1 

3  | 

EMC: 

IvO 

1 

2  1 

EMC; 

1 

1 

/•  too  HEtSOSEHEST  EEkOES  •/ 

1V1 

1 

1  | 

DO  1-0  TO  Hi 

1v  2 

1 

2  1 

X E S  (It  II  -XEE  (1,1)  •H,*2i 

IV  1 

1 

2  1 

E  H  C  # 

1 

1 

/•  COMPUTE  ITXES  V 

1»  i 

1 

1  | 

HTXE  E-0 ; 

1V5 

1 

1  | 

00  1-5  TO  H; 

1V6 

1 

2  1 

00  J-G  TO  Hi 

1*7 

1 

3  1 

STXES  (I)  -STXEE(I)  ‘EE  (J|  ‘XEE(J.I)  i 

ivg 

1 

3  1 

ehc: 

1V9 

1 

2  1 

END; 

1 

1 

/•  COMPUTE  Oil  •/ 

150 

1 

1  l 

Dll *0; 

161 

1 

1  1 

00  I-G  *0  Hi 

152 

1 

2  1 

311'Cll‘HTXEE  II)  ‘EE  (I)  i 

15  1 

1 

2  1 

EHC; 

15v 

1 

1  | 

012*8TI‘SQPTIDtt)  /STS; 

155 

1 

1  | 

IP  11>0  THEM  GO  TO  DD; 

156 

1 

1  | 

IP  |StX(G)>90)6{Stl(G)  <270)  THEE  SWj  ELSE  11*1; 

158 

1 

1  | 

IP  Sll  fG)  <180  THIS  12*1;  ELSE  12>2i 

160 

1 

1  | 

IP  ESI  Z  (HI  >90)1(511(8)  <270)  T<EI  13-2;  ELSE  13*1; 

162 

1 

1  | 

IP  51!(H)<1S0  THEE  11-1;  ELSE  1«*2| 

16V 

1 

1  • 

ELC1-0 i 

165 

1 

1  , 

DO  I-G  TO  H; 

166 

1 

2  1 

IP  El(I)>ELCl  THEE  ELC1«EL(I)' 
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n/1  OPTINIIING  CONPILEN  AIAL:  PNOCIDUNI  OPTIONS  (mill)  ; 


ST*? 

lev 

nt 

U7 

i 

2 

END; 

( 

16H 

i 

1 

Pin  5«  :p(2)  ; 

l 

1 

| 

169 

i 

1 

P'JT  EDIT  (N  |A1)  .E  (A2)  ,*  TO  '  ,  N  (A 3)  ,  E|ll|  ,  •  PASS  ELEVATION  AT  CA  IS* 

1 

1 

1 

,ELC»,  •  OEG.  ')(l  (101  ,({!).>  15,1)  ,A) ; 

) 

170 

i 

1 

put  s«:p(2i  i 

1 

171 

1 

1 

PUT  EPITI'ANT*, ‘ERROR  IN  SECONDS  PEN  VARIABLE*  (  *  RBS  *| 

\ 

} 

(1  (  111  ,»,I  (9)  ,A,I  (6f  (  ,A(  ; 

( 

1 

172 

i 

1 

PUT  SN  IP  EDIT  <*Al*,*LAT*,*LON*  ,*HT*,*IS*,*  IS*  ,  • IS* , *D • , *DH* . • . 

( 

•TOTAL*! 

( 

l  ) 

17] 

i 

1 

(I I  111.  A,  I  (6)  ,2  (A.  I  (•)),«  (A,  MR)),  A,  I  (10)  .A,  1(9)  ,A,1(S),A>< 

PUT  SKIP; 

( 

17* 

i 

1 

DD:  PUT  SNIP  EDIT  (AI,NTS»ABS(«TN»rN/»TI)  ,»TS»R/SmT(»T»)  .DAI) 

) 

( 

.*  | 

(I(11(.P(2l.S(T).T0(I<J)  ,■  ; 

) 

1 7  S 

i 

1 

END;  /MND  AIV 

( 

176 

i 

0 

PUT  SNIP  («<  J 

177 

i 

0 

GO  TO  NET; 

1 

1 

178 

i 

0 

END  ARAL; 

1 

) 
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